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Abstract 


Our  primary  aim  is  to  "^build^versions  of  generalised  Gaussian  pro¬ 
cesses  from  simple,  elementary  components  in  such  a  way  that  as  many 
as  possible  of  the  esoteric  properties  of  these  elusive  objects  become  in¬ 
tuitive.  For  generalised  Gaussian  processes,  or  fields,  indexed  by  smooth 

l\  <u  r  -i 

functions  or  measures  ouiK  our  building  blocks  will  be  simple  Markov 

* 

processes  whose  state  space  is  .  Roughly  speaking,  by  summing  func¬ 
tions  of  the  local  times  of  the  Markov  processes  we  shall,  via  a  central 
limit  theorem  type  of  result,  obtain  the  Gaussian  field. 

This  central  limit  result,  together  with  related  results  indicating 
how  additive  functionals  of  the  Markov  processes  generate  additive  func¬ 
tionals  of  the  fields,  yield  considerable  insight  into  properties  of  gener- 
alised  Gaussian  processes  such  as  Markovianess,  self-similarity,  "locality 
of  functionals,  etc. 

Although  the  paper  is  comprised  primarily  of  new  results,  and  de¬ 
spite  the  fact  that  the  subject  matter  is  somewhat  esoteric,  our  aims 
are  primarily  didactic  and  expository  we  want  to  try  to  initiate  the 
uninitiated  into  some  of  the  mysteries  of  generalised  processes  via  an 
easily  understood  model. 
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1.  Introduction. 

This  paper  is  about  generalised  Gaussian  processes  and  «onie  of  their  properties,  or, 
to  be  more  precise,  it  is  about  a  way  to  think  about  these  things  without  becoming  too 
confused.  Generalised  Gaussian  processes  can  be  succinctly  defined  as  isometric  mappings 
from  a  Sobolev  space  to  a  space  of  Gaussian  variables,  and  virtually  all  their  interesting 
properties  can  be  stated  in  terms  of  the  associated  Fock  space.  While  such  definitions  and 
statements  are  both  neat  and  precise  they  do  nothing  to  help  the  non-expert  understand 
what  these  rather  esoteric  objects  are.  Our  plan  is  to  write  for  the  novice  -  we  shall  start 
at  a  point  which  seems  to  us  to  be  a  natural  beginning,  and  follow  a  natural  sequence  of 
observations  and  constructions  that  will  ultimately  show  Sobolev  and  Fock  spaces,  along 
with  much  that  goes  with  them,  in  a  fashion  that  is  intuitive  for  anyone  who  knows  what 
a  Markov  process  and  a  central  limit  theorem  (CLT)  are. 

Since  generalised  Gaussian  processes  and  all  that  goes  with  them  are  actually  a  part 
of  Mathematical  Physics,  (they  are  “just"  Euclidean  quantum  fields,  but  we  shall  take 
up  this  point  later)  and.  as  such,  are  supposed  to  describe  the  behaviour  of  elementary 
particles,  let  us  start  by  watching  such  a  particle.  In  fact,  let  us  take  .V  >  0  such  particles, 
each  one  of  which,  at  some  common  and  fixed  point  of  time,  "pops'’  into  existence  and 
then  proceeds  to  wander  about  space,  which  we  take  to  be  R'1,  d  >  2.  Assume  the 
"birthplaces”  of  these  JV  particles,  xi,...,z.v  are  distributed  randomly  in  R'1  according 
to  some  law,  and  that  each  particle  executes  a  Markov  motion  through  R'*,  described  by 
A'  independent  Markov  processes  A'i(0-  •  ■  • ,  A'v(f).  t  >  0.  To  make  our  story  sound  more 
like  Physics,  we  shall  also  assume  that  each  of  these  particles  begins  its  life  with  a  charge 
<r.,  where  P(<r,  =  +1)  =  P(<x,  =  -1)  =  1/2.  The  charges  of  different  particles  are  assumed 
independent. 

To  complete  the  story,  suppose  that  at  time  t  the  original  unit  charge  has  dissipated  by 
a  factor  of  e_(,  and  that  as  each  particles  performs  its  Markov  motion  through  RJ  it  leaves 
each  point  it  hits  with  a  charge  equal  to  its  own  charge  at  the  time  of  hitting  (i.e. 

What  we,  as  amateur  physicists,  would  like  to  know  is  what  does  the  charge  distribution 
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throughout  RJ  look  like  after  all  our  particles  have  decayed  to  inconsequentially.  Clearly, 
as  prob  ibilists.  we  shall  only  be  interested  in  the  .V  — *  oo  limit  of  this  question! 

Considering  only  one  such  particle  for  the  moment,  it  is  clear  that  the  key  to  solving 
our  problem  must  lie  in  the  definition  of  some  sort  of  weighted  local  time ,  a  random  function 
LZ(.Y)  that  would  measure  the  amount  of  time  the  path  X(t).t  >  0.  spends  at  the  point 
x  €  RJ.  weighted  by  its  charge  when  visiting.  That  is,  we  would  like  to  make  sense  of  the 


equation 


/•OO 

Lx{X)  =  I  e~t6(X(t)- x)dt, 
Jo 


where  6  is  the  Dirac  delta  function.  For  the  moment  let  us  take  (1.1)  at  face  value,  even 
though,  as  we  shall  see  later,  it  has  no  fully  rigourous  meaning.  Then  our  problem  reduces 
to  studying  the  behaviour  of  the  sum 

,v 

Q,\(x)  :=  N~lf2  ^2ffkLx(Xk) 

k  —  l 

as  .V  —  oc.  If  there  are  enough  moments  about  (and  we  shall  ensure  that  there  are)  then 
the  independence  of  the  A*  should,  with  a  standard  CLT,  ensure  that  the  9y  converge  to 
a  Gaussian  random  field  on  Rd.  This  would  be  a  first  step  towards  solving  our  problem. 
(But  only  a  first  step,  since,  as  amateur  physicists,  we  would  also  like  to  know  something 
about  how  particles  interact.  So  far,  nothing  in  our  model  allows  for  interaction.) 

Unfortunately,  this  simple  minded  approach  has  one  basic  flaw  in  it  -  the  integral  in 
(1.1)  does  not  generally  exist.  Probabilistically  this  is  because  most  Markov  processes  (e.g. 
Brownian  motion)  do  not  hit  points  in  R d  d  >  2,  with  non-zero  probability.  (In  terms  of 
the  physics  of  the  limit  random  field  this  turns  out  to  be  related  to  such  elegant  concepts  as 
the  Heisenberg  uncertainty  principle!)  To  get  around  this  difficulty,  we  note  that  instead 
of  asking  how  much  weighted  time  each  particle  spends  at  a  given  point,  we  could  ask 
about  the  time  it  spends  in  some  small  set  A.  Writing  1.4  for  the  indicator  function  of  A, 
this  would  lead  us  to  a  weighted  occupation  measure  which  we  could  denote  by 


La(X)  =  L[1a:X) 


-  I 


e-‘l  A(X(t))dt. 
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There  is  no  question  that  this  object  exists  without  anything  but  the  mildest  measurability 
conditions  on  A”.  Of  course,  if  (1.1)  were  justifiable,  then  we  would  have 


Z(1a:A)-  [  LA*)djr=  [  lA(x)Lx(X)dx. 

Ja  J  r“ 


Freely  generalising  these  equivalences,  let  us  now  take  /  to  be  any  “nice"  (we  shall  be 
more  specific  later)  function  on  R1*,  and  consider  a  time  and  space  weighted  local  time  of 
the  form 


(1-3) 


L(/;X)  :=  /  f(x)Lx(X)dx. 

Jnd 


With  /  an  indicator  function  we  recover  the  well-defined  (1.2).  With  /  a  delta  function  we 
recover  the  ill-defined  (1.1).  With  /  belonging  to  just  the  right  class  of  functions  (Schwartz 
or  Sobolev  space  -  but  wait  for  details)  (1.3)  can  be  shown  to  make  perfect  sense.  Another 
way  to  write  (1.3)  is  to  notice  that  if,  in  fact,  Lx  existed,  then  we  would  have 


(1.4) 


L(f 


;-Y)=  / 

Jo 


e-7(-V(0)  dl- 


In  this  formulation  L[f:  X)  is  certainly  well-defined  with  only  minimal  assumptions  on  /. 
Nevertheless,  (1.4)  is  undesirable  for  two  reasons.  Firstly,  it  does  not  generalise  easily  to 
the  measure-indexed  local  times  we  shall  consider  in  a  moment,  and,  secondly,  it  does  not 
have  the  intuitive  appeal  of  (1.3). 

If,  in  (1.3),  we  think  of  /  as  being  a  weighting  function  for  the  point  indexed  local 
time,  then  it  is  only  a  small  step  to  think  of  weighting  by  measures.  That  is.  if  n  is  a  nice 
measure  on  RJ  then  we  can  make  sense  of  the  weighted  local  time  L(p.  A)  which  we  write, 
symbolically,  as 

I(/i;-Y):=  /  LAX)n(dx). 

JR* 

Again,  although  Lx  may  not  exist,  it  is  possible  to  make  sense  out  of  I(^.A'). 

Returning  now  to  our  CLT.  it  is  clear  that  the  sums  we  should  really  be  working  with 
are  either 


15) 


*v(/)  :=  -V-,/2I>£(/..Y< 

k~l 


■  ^  -  f  -  «  ^1.  ■  -  <  .S  1.1. 


t  —  ■*-  ■*—  _  i  *  -  *  *  * 


•y  > 


(1.6) 


i.  =  l 

The  limits  of  these  sums  will  he  Gaussian  variables  parameterised,  as  we  vary  the  weights 
/  and  /<.  by  either  a  family  of  functions  or  a  family  of  measures.  The  Gaussian  processes 
thus  obtained  will  be  the  generalised  Gaussian  processes  that  we  mentioned  in  our  opening 
sentence,  and  the  CLT  described  will  be  our  tool  for  studying  them.  This,  however,  does 
uot  represent  the  end  of  our  task,  for  we  have  not  yet  found  a  way  to  describe  interactions 
between  particles. 

For  the  sake  of  simplicity,  let  us  consider  for  the  moment  only  a  very  simple  type  of 
interaction.  Recall  that  as  each  of  our  particles  passes  a  point  x  6  RJ  it  leaves  behind  a 
charge  ae~t .  where  t  is  the  time  at  which  it  hits  x.  Suppose  another  particle  hits  x  at  time 
s.  leaving  behind  its  charge,  and  that  charge  behaves  in  a  multiplicative  fashion  (a  highly 
unrealistic  assumption,  but  we  are  only  amateur  physicists).  That  is,  the  charge  left  at  x 
from  the  interaction  of  A',  and  Xj  hitting  at  times  t ,  and  t}.  respectively,  is 
Consequently,  if  we  are  interested  in  the  total  interaction  charge  at  x  then  we  must  study 
a  new  intersection  local  time  of  the  form 


e  -  x)e  *8(X,(8)-x)dtd8. 

Of  course,  if  (1.1)  with  its  single  delta  function  doesn't  exist,  then  a  fortiori,  neither  will 
this  new  intersection  local  time.  Nevertheless,  we  can  proceed  much  as  before,  firstly 
replacing  the  delta  function  in  (1.7)  with  the  indicator  function  of  some  set,  and  then  with 
a  general  function  or  measure,  to  ultimately  suggest  a  weighted  intersection  local  time 
which  could  be  symbolically  written  as 

(1.8)  L(/i;  A',..Y;)  :=  /  Lj(Xt,Xj  )n(dx). 

Jr* 


1.7)  Lx(X,,X,):=  i  [ 

Jo  Jo 


Working  somewhat  harder  than  we  have  in  the  above  argument,  we  shall  see  later  that  the 
left  hand  side  of  (1.8)  can  be  given  a  rigourous  meaning. 


With  a  notion  of  intersection  strength  (between  two  particles  at  a  time)  now  detiued. 
we  could  hope  to  study  a  total  intersection  process  by  looking  at  the  sum 

(19)  *.v(p)  :=  A',). 

the  sum  being  over  all  i  <  j,  i.j  =  1 . .V.  As  probabilists  we  would  now  like  to  send 

-V  — •  oc  in  this  expression.  It  is  clear  that  no  simple  CLT  will  work  here,  since  the 
summands  in  ( 1.9)  are  dependent.  Nevertheless,  one  could  hope  that  since  looking  at  the 
local  times  of  the  Markov  paths  themselves  lead  us  to  a  generalised  Gaussian  process,  then 
looking  at  a  functional  of  these  paths  (their  intersections)  should  lead  us  to  a  functional 
of  the  Gaussian  limit.  In  fact,  this  is  precisely  what  happens,  and  (1.9)  converges  to  what 
is  known  as  the  Wick  square  of  the  basic  Gaussian  limit.  The  Wick  square  lives  in  the 
first  interesting  part  of  Fock  space,  and  so  we  shall  now  have  a  way  of  thinking  about 
generalised  Gaussian  processes  and  their  Fock  spaces  in  terms  of  simple  Markov  particles. 
This  is  precisely  what  we  are  after! 

Before  we  can  go  any  further,  we  are  going  to  have  to  start  making  the  above  picture 
somewhat  more  precise  by  defining  our  terms  more  carefully.  In  the  following  section  we 
shall  therefore  collect  some  information  on  Markov  processes  and  their  additive  function¬ 
als.  (Both  the  local  times  L(p; A',)  and  L(p;  A', , X})  are  additive  functionals.)  In  Section 
3  we  shall  do  the  same  for  generalised  Gaussian  processes  and  their  additive  functionals, 
(e  g.  Wick  squares)  In  Section  4  we  shall  then  be  able  to  continue  where  the  introduction 
leaves  off,  by  properly  formulating  general  limit  theorems  that  allow  us  to  describe  Gaus¬ 
sian  processes  and  their  additive  functionals  via  Markov  processes  and  their  functionals. 
Following  this,  we  consider  the  insights  that  these  results  lead  to  in  Sections  5  and  6. 

The  casual  reader  can  then  stop,  for  the  harder  work  starts  in  Section  7  with  a  proof 
of  special  cases  of  the  main  result  via  the  method  of  moments.  (This  is  not  too  hard.)  In 
order  to  prove  the  most  general  form  of  the  main  result,  we  need  to  know  a  little  more 
about  the  multiple  Ito-W'iener  integrals  that  describe  Fock  space  than  what  is  given  in 
Section  3.  Section  8  provides  this,  and  also  contains  the  final  proof,  which  is  also  not  too 


hard,  given  that  everything  has  been  carefully  set  up  for  it.  Some  closing  comments  are 
made  in  the  concluding  Section  9. 

Before  commencing,  however,  we  have  some  general  historic  comments  and  some  ac¬ 
knowledgements  to  make.  Our  program  of  viewing  Gaussian  fields  from  a  Markov  outlook 
is  not  new.  The  fields  that  we  shall  consider  are  known  to  mathematical  physicists  as 
Euclidean  quantum  fields,  and  are  imaginary  time  versions  of  real  quantum  fields.  As  such 
they  have  been  studied  extensively.  Furthermore,  it  has  long  been  known  that  they  are 
related  to  Markov  processes,  the  seminal  work  in  this  area  being  due  to  Symanzik  (1969). 
This  approach  has  been  taken  up  and  developed  at  length  in  a  program  of  Dynkin’s  (1980, 
1983,  1984a, b,  for  example)  to  relate  Gaussian  and  Markov  processes.  However,  through¬ 
out  all  these  projects  there  is  no  attempt  to  tie  in  the  Gaussian  and  Markov  processes  via 
a  physical  model  of  the  kind  described  above.  The  only  exception  to  this  is  in  the  work  of 
Wolpert  (1978a.b),  which  builds  a  model  very  similar  to  that  we  have  constructed.  In  fact, 
it  was  Wolpert's  elegant  construction  that  lead  us  to  the  results  that  follow.  Our  results 
differ  from  his  in  two  ways.  Firstly,  they  cover  a  much  wider  field  of  models,  (he  deals 
only  with  the  so-called  "free  field"  and  Brownian  motion  as  the  Markov  process)  and  we 
deal  in  considerable  more  depth  with  the  consequences  of  the  model.  Secondly,  our  style 
of  proof  is  generally  quite  different,  and,  we  feel  much  neater. 

The  main  contribution  of  our  paper,  then,  is  not  so  much  the  introduction  of  a  com¬ 
pletely  new  idea  b>’*  rather,  the  development  of  a  way  of  looking  at  things  and  an  expo¬ 
sition  of  that  way  that  we  hope  will  be  at  a  level  to  make  it  available  to  a  wide  audience. 

Regarding  our  audience  -  since  a  substantial  part  of  our  aim  has  been  didactic  we 
have  not  written  this  paper  for  the  expert  quantum  field  theorist  who  already  has  a  well- 
developed  intuition  about  generalised  Gaussian  processes.  (Although  we  do  hope  that  even 
he  may  find  something  of  interest  here.)  Rather  we  have  written  with  the  curious  novice 
in  nnnd.  In  particular,  we  should  point  out  that  most  readers  will  find  either  Section  2 
or  Section  3  “well  known",  although  we  expect  that  few  will  feel  that  way  about  both. 
Feel  free  to  skip  the  section  that  you  already  know.  The  paper  also  turned  out  much 


longer  than  we  had  ever  expected  it  would.  After  all,  our  avowed  aim  has  been  to  produce 
something  readable,  and  not  a  forbidding  treatise.  However,  after  a  little  practice  the 
reader  will  learn  to  separate  the  interesting  and  important  parts  from  the  clumsy  and 
space  consuming  notations  that  seem  to  be  endemic  to  this  area  of  Probability,  and  the 
reading  of  it  all  should  go  much  faster  than  one  would  expect  in  the  beginning. 

Finally,  the  acknowledgements.  Our  dept  to  Professor  Dynkin's  recent  work  will  be 
obvious  by  the  citations  to  it  throughout  the  paper.  What  is  not  obvious  is  the  debt  we 
owe  him  for  an  inspiring  set  of  lectures  that  he  gave  during  a  visit  to  Israel  in  1983  which 
did  much  to  motivate  our  work.  It  is  a  pleasure  to  acknowledge  that  debt  now.  Haya  Kaspi 
spent  a  substantial  amount  of  time  answering  a  lot  of  questions  about  Markov  processes 
in  the  early  stages  of  this  work.  We  could  not  have  managed  without  her.  Murad  Taqqu 
helped  us  to  understand  parts  of  Dynkin  and  Mandelbaum  (1983)  and  Mandelbaura  and 
Taqqu  (1984).  Most  of  this  paper  was  written  while  both  authors  were  visiting  the  Center 
for  Stochastic  Processes  in  Chapel  Hill.  The  hospitality  of  the  Center  and  its  directors 
was.  as  always,  gracious  and  most  appreciated. 


2.  Markov  Processes  and  their  Functionals. 

2. a.  Markov  processes  and  systems. 

The  treatment  of  Markov  processes  and  their  functionals  given  in  this  section  comes, 
essentially,  from  Dynkin  (1981),  to  which  we  refer  the  reader  for  missing  detail. 

Let  X(t).t  >  0,  be  a  Markov  process  on  a  probability  space  (0,7,P),  taking  values 
in  Rd.rf  >  1.  We  assume  the  existence  of  a  stationary,  symmetric  transition  density 


(2.1] 


pt(x,y)  =  pt(y,  *) 


satisfying  JRd  Pi{x,y)dy  =  1  for  all  x.  For  each  “initial  point”  x  the  “probability  starting 
at  /"  is  given  by 


. XilteBk}=  /  ...  PtAx'dyl)pti-tAVi'dy2)..  Ptk-tt-AVk-iidyk) 

J  B  i  JBk 


while  each  <r-fmite  initial  measure  m  defines 


(2.2) 

For  0  >  0  the  function 


P-(.);=  / 

JR‘ 


Px  (•)m(dx). 


(2.3) 


rOO 

Q9(x.y)  ■=  /  e~0tPt(i, 
Jo 


y)dt. 


is  called  the  Green  function  for  A'  or  pt .  To  obtain  a  heuristic  feel  for  the  Green  function, 


note  that  for  an  arbitrary  set  A  C  Rd 


/  ge(*,v)l-A(v)dy=  [  [  e  etPt(x.y)iA(y)dydt 

J  It'  J  Rd  Jo 


(2.4) 


e~9tpx{X(  € 


=  E1  J 
Jo 


e~etlA(Xt)dt. 


But  this,  with  0  =  1,  is  the  expected  value,  when  A  starts  at  the  point  r,  of  the  weighted 
occupation  measure  of  (1.2).  Replacing  lA  by  a  delta  function  shows  that  g9(r.y)  can 


be  interpreted  as  the  expectation  of  the  exponentially  weighted  local  time  of  (1.1).  with 
weighting  e~6t  rather  than  e~* . 

In  what  follows  we  shall  generally  need  to  consider  the  Markov  system  of  k  >  1 


independent  copies  A'1 . Xk  of  X.  which  we  write  as  a  Revalued  function  on  R^  via 

(2.5)  X|t)  :=  ( -V 1  ( / 1 ) . Xk(tk)). 


For  notational  convenience,  we  introduce  fc-dimensional  “transition"  and  "Green"  functions 
as 

k  k 

Pt(x.y)  :=  Yl  Pr{*'-y')-  /(x.y)  :=  II 

t=i  t=i 

with  x  =  (/l, _ xk)  6  RJfc.  Note  that  by  (2.1)  we  always  have  ge(x.  y)  =  gd(y,x). 

Throughout  the  remainder  of  the  paper  we  shall  assume  that  all  processes  are  right 
in  the  sense  of  Dynkin  (1981). 

2.b.  Additive  functionals. 

To  start  with  an  easy  case,  let  A'  be  a  Markov  process  in  RJ  and  X  a  system  of  k 
independent  copies  of  it.  Let  6  be  a  positive,  bounded,  measurable  function  on  RJfc  and 
0  >  0.  Then  the  path  integral 

(2.6)  F6(I):=  + ■+tk]b(X{t))dt 

defines  a  family  of  random  variables,  indexed  by  open  intervals  I  C  R^_,  with  the  following 
properties: 

(2.7)  For  every  path,  the  function  Fe  can  be  extended  to  a  measure  on  R^, 

(2.8)  For  every  I,  Fe (I)  is  an  element  of  the  minimal  <r-algebra  generated  by  the  X(t).  t  €  I. 

(2.9)  A  shift  t  —  s  +  t  of  the  path  of  X  induces  an  analogous  shift  of  the  corresponding 
measure. 

The  fact  that  we  have  written  the  exponential  term  in  (2.6)  explicitly,  and  not  as  part 
of  the  function  b.  is  to  recall  the  exponential  weakening  of  the  Introduction.  However,  in 


9  M.iPJP.'PJP  J  PJ » J  »J  •  Y  .■  •  ■  « « UW'TW'1  ^  ? 


T*T 


one  form  or  another,  some  damping  is  necessary  for  the  integral  to  be  finite.  (Even  the 
simple  case  b  =  1 A  would  lead  to  a  divergent  integral  if  A'  is  neighbourhood  recurrent, 
unless  the  exponential  is  included.)  Rather  than  damping  X.  we  could  leave  it  untampered 
until  some  exponential  killing  time.  This  would  lead  to  a  similar  theory  and.  in  terms  of  the 
Introduction,  w’ould  take  us  closer  to  Wolpert’s  model.  Purely  for  reasons  of  mathematical 
convenience,  we  prefer  to  dampen  in  a  continuous  fashion. 

Any  family  with  the  properties  (2. 7) -(2. 9)  is  called  a  homogeneous  additive  functional 
of  X.  The  class  of  all  such  functionals  is  much  larger  than  that  obtained  by  integrals  of 
the  form  (2.6).  In  fact,  what  we  shall  be  most  interested  in  will  be  situations  in  which 

b(X(t))“  ="  <l>fc-l(A'l(t)  -  A'2(t) . Xk~l  (t)  -  A'fc(t)).  with  fim  the  Dirac  delta  function 

on  F"1"*.  In  such  a  case,  the  additive  functional  Fe  measures  the  set  of  common  points  of 

A'1 . Xk .  in  a  sense  to  become  clearer  later. 

To  extend  (2.6).  let  7.  v  be  measures  on  RJlc,  and  introduce  the  inner  products 

=  /  -,(dx)/(x.y)i/(dy) 

(2.10)  J 

=  J  -;l dxl . . dyk). 

Let  Mfl  k  =  Mk(gt> )  denote  the  set  of  all  <r-fini*e  7  for  which  (7,7)0  <  00.  Also,  for 
B.t  €  Rk  write  g  <  t  if  s'  <  t'  for  all  i,  and  if  s  <  t  then  let  (s.t)  denote  the  open 

interval{u  .  s‘  <  u'  <  t'.i  —  1 . Is}.  Then  the  following  theorem,  linking  measures  in 

A!*1  -k  and  additive  functionals,  is  a  trivial  extension  of  Dynkin  (1981).  (Dynkin  treats  the 
case  0  =  1.  The  extension  to  general  0  is  immediate.)  For  notational  convenience,  we  write 
|t|  for  <!+...  +  and  ||t||2  for  I2  +  ...  +  <£. 

Theorem  2.1.  Let  the  k  components  of  X  be  independent  copies  of  a  symmetric  right 
process  with  initial  measure  rn  =  Lebesgue  measure.  Then  to  every  measure  7  6 
there  corresponds  an  additive  functional  Ft  of  X  with  the  following  properties: 

I  There  exists  a  negligible  set  fl0  such  that  for  all  w  Ho  F!j(u’;  ( s.  t ) )  is  finite  for  all 
0  <  s  <  t  and  is  continuous  in  s  and  t 
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II:  For  every  positive,  Bore 1  f  :  R*j_  x  Rdk  — *  ft 


(2.1i)  Em  f{t,X{t))F*(dt)  =  /  e-e^/(t,yh(dy)dt 

J R*  J H*  JRik 

Whereas  Theorem  2.1  guarantees  us  an  F 9  for  each  7  6  M9'k,  one  can  also  go  the 
other  way.  That  is,  given  an  additive  functional  F  of  X,  the  so-called  spectral  measure  7 f 
defined  by 

7 F(A)  -.=  0Em  j  lA(X(t))F(dt) 

(A  C  Rdk)  defines  a  measure  satisfying  (2.11). 

For  a  measure  7  €  Me'k  the  corresponding  additive  functional  F9  can  be  constructed 
in  one  of  two  ways.  If  7  is  absolutely  continuous  with  respect  to  Lebesgue  measure  on 
Rdk,  so  that  7(dx)  =  6(x)dx  with  J  6(x)j/9(x,  y  )6(y)dxdy  <  00,  then  F9  can  be  defined  as 

(2.12)  F9(/)  =  J^e~9^b(X(t))dt,  I  C  R*, 

so  that  we  return  to  the  form  of  the  simple  functional  of  (2.6).  In  the  general  case.  F9 
is  defined  as  a  limit  of  path  integrals  of  this  form,  with  the  function  b  being  densities  of 

smoothed  versions  of  7.  To  be  more  precise,  for  0  :=  (0i . ),  let  such  a  density  be 

given  by 

(2.13)  b9s{x)=[  e~eltlpt(x,y)-f{dy) 

Jr1"1 

so  that  the  smoothed  versions  of  7  given  by 

T  g(dx)  =  b9s(x)dx 

satisfy  (7  -  7^.7  -  7^)5  —  0  as  J|0j|  — >  0.  Then  if  7  is  finite  (i.e.  7(Rdfc)  <  oc)  the 

F’f  (/):=  j 


functionals 


converge,  in  C2{P ),  to  the  functional  F '*  of  the  theorem.  If  7  is  not  finite,  but  belongs  to 
Md,k.  then  there  exists  a  finite  measure  -  '  and  a  function  h  such  that 


(2.15)  7 (dx)  =  h{x)~j'(dx ),  7'  €  Xs,*\ 

In  this  case,  F *  is  defined  as  F*(dt)  =  /i(X(t))Fy(dt). 

2.c.  An  example. 

To  gain  some  feeling  for  the  above  results,  and  as  a  precursor  to  Section  5,  we  consider 
one  example.  Take  k  =  d  =  2,  6  =  1,  and  A'i  and  A \  independent  Brownian  motions  in 
R2.  Then  for  x.  y  G  R2 

gl{x.y)=f  e-,(2jr/)-1  exp{-~]|x  -  v\\2 /0 

(2.16)  Jo 

=  eA'0(i|/-y((), 

where  A'o  is  a  modified  Bessel  function  and  c  an  uninteresting  constant.  For  A  C  R2,  take 
7.4  to  be  Lebesgue  measure  on  the  two-dimensional  set  DA  in  R4  given  by 

(2.17)  Da  =  {(/i./2./3,/4)  :  (/1./2)  €  .4,(/j,/4)  €  A,/!  =  *,.*2  =  /<}. 

Then,  since  A'0(c)  decays  exponentially  fast  for  large  z,  and  has  only  a  logarithmic  sin¬ 
gularity  at  :  =  0,  it  is  easy  to  check  that  (7.4,  7.4)1  <  oc  and  so  7.4  €  X12.  Thus  an 
F71a  satisfying  (2.11)  exists,  and  from  the  form  of  7.4  and  (2.11)  it  is  relatively  easy  to  see 
that  F*a  charges  only  those  t  =  ( f j ,  f2 )  for  which  A'tJfj)  =  A'2(f2)  and  A't(<,)  6  A.  The 
integral  f  F*A(d t)  now  gives  us  a  candidate  for  the  weighted  intersection  local  time  of  A'i 
and  A'2  while  they  are  in  A.  Replacing  A  by  R2  and  7.4  by  a  measure  n  living  on 
and  satisfying  (/i,ji)i  <  00  yields 

L(iv.  A'.,A'a)=  f  F^(dt) 

J  R* 

as  a  properly  defined  version  of  (1.8). 

Thus,  we  now  have  the  basic  components  for  the  limit  theorem  discussed  in  the  In¬ 
troduction,  and  so  we  turn  to  a  discussion  of  what  will  ultimately  yield  the  limit  process. 


S.  Gaussian  Fields  and  their  Functionals. 

S.a.  Gaussian  random  fields. 

Let  g(i.y)  be  a  positive  definite  function  on  Rd  x  Rd.  Then,  in  its  simplest  form, 
the  Gaussian  random  field  with  covariance  function  g  is  the  family  of  Gaussian  random 
variables  *(/),/  €  RJ  with  zero  mean  and  covariances 

(3.1)  £{$(*)$(*/)}  =  g(x<y). 

(Note  the  sequence  of  ideas  -  for  us  the  covariance  function  comes  first  and  the  random 
field  second.)  The  random  fields  that  arise  in  quantum  field  theory,  and  the  interesting 
ones  that  arise  as  the  limit  processes  mentioned  in  the  Introduction,  correspond  to  positive 
definite  g  with  g(x, /)  =  oc  for  all  x  €  Rd.  giving  the  point  indexed  field  $(x)  of  (3.1) 
infinite  variance  and  so  leaving  it  poorly  defined.  The  usual  way  around  this  is  to  change 
the  parameter  space  of  $  to  a  space  of  functions.  There  are  two  natural  choices.  One  is  Sj, 
the  Schwartz  space  of  C ^  functions  that  decrease  at  infinity  faster  than  any  polynomial. 
Given  a  rich  enough  probability  space  (f),  /,  P)  the  Gaussian  field  on  Sj  with  covariance 
kernel  g  is  defined  as  the  continuous  linear  mapping  from  Sd  to  the  zero  mean  Gaussian 
variables  in  £°(P),  (with  the  topology  on  the  latter  space  being  that  of  convergence  in 
probability),  such  that  the  covariance  between  two  values  of  the  field  is  given  via  the 
bilinear  functional 

(3.2)  B9(f,h)  =  £{*(/)*{/.»  =  j  I  f(x)g(x,y)h(y)dxdy. 

(«  S!{  is  the  topological  dual  of  Sd  -  i.e.  the  space  of  generalised  functions  of  temper¬ 
ate  growth,  or,  tempered  distributions  -  $(•)  can  be  considered  as  a  Sj-valued  random 
variable.  Hence  the  now  almost  archaic  term  '‘generalised  Gaussian  process"  used  in  the 
Introduction,  a  term  which  we  now  drop  in  favour  of  "random  field".) 

Clearly,  if  (3.2)  is  to  make  sense,  we  require  Bg(f. /)  <  oc  for  all  /  €  S,t.  so  we  are 
limited  in  our  choice  of  covariance  kernels  g  to  those  for  which  this  is  true.  In  order  to 
allow  the  definition  of  fields  for  arbitrary  g.  we  restrict  the  parameter  space  to  the  Sobolev 
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space  S(g)  of  C ^  functions  of  finite  norm,  where  the  inner  product  in  S(<7>  is  given  by 
=  J  J  f(f)g(x.y)h(y)didy.  A  Gaussian  field  on  S(t?)  is  defined  as  for  one  on  Sj. 
with  covariance  functional  (3.2). 

(Note  that  here  Bg(f.h )  =  ( f.h)g ,  leading  to  our  opening  comments  about  Gaussian 
fields  being  isometric  mappings.) 

If  a  point  indexed  field  satisfying  (3.1)  exists,  then  we  can  always  create  a  5j  or  S ( j?) 


indexed  version  of  it  by  setting 


*(/)  =  iumx)dx. 


and  (3.2)  will  be  satisfied.  However,  function  indexed  fields  satisfying  (3.2)  often  exist, 
even  when  they  cannot  be  represented  in  the  above  form. 

Another  way  to  index  random  fields  is  by  the  family  of  measures  introduced  in  the 
previous  section,  i.e.  Ml  =  Ml(ff)  =  {v  '■  (fb/i }g  <  oc},  where  we  had 


J  J  H(dx)g(x.y)u{dy), 


for  any  two  Borel  measures  (t.v.  The  corresponding  Gaussian  field  maps  MMff)  to  zero 
mean  Gaussian  variables  with  £  {$(/t)$(t/) }  =  (y,v)g.  As  above,  if  $(/)  exists,  we  have 
$(/i)  =  J  ${i)n{djr).  This  index  set  seems  to  have  been  initially  considered,  in  the  context 
of  the  problems  that  will  interest  us,  by  Albeverio  and  Hoegh-Krohn  (1979),  and  was  then 
taken  up  by  Dynkin  (1980).  Since  it  turns  out  to  be  the  most  natural  parameter  space 
over  which  to  study  the  Markov  property  of  random  fields  we  shall  concentrate  on  it  in  the 
future.  However  it  is  important  to  note  that  since  measures  in  Ml(g)  can  be  approximated 
by  elements  of  Sd  or  S(p),  and  vice  versa,  the  theories  associated  with  all  three  parameter 
spaces  are  essentially  the  same,  (except  for  their  Markov  properties  -  for  details  on  this, 
see  Albeverio  and  Hoegh-Krohn  (1984)  and  references  therein). 

One  example  of  particular  interest  is  the  so-called  free  field  of  mass  m  >  0,  for  which 


the  covariance  kernel  is  given  by 


gm[x.y )  —  3m  \x 


y]  (2x)d  L 


t  m2  +  j|fc||2 
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Since  g(i.  .r)  =  x  for  d  >  2.  the  free  field  does  not  exist  pointwise,  hut  since  the  singular¬ 
ities  of  g  are  then  at  worst  polynomial,  free  fields  on  S,i.  S(f7)  or  Ml(tf)  can  be  defined. 

The  kernel  (3.5)  can  be  looked  upon  in  a  number  of  ways.  It  is  the  kernel  of  the 
operator  (-A  +  m2)-1  in  Rd,  afact  important  in  Euclidean  field  theory.  More  importantly 
for  us.  however,  is  that  it  is  also  the  Green  function  g 6  of  Brownian  motion  in  Rd  with  0  = 
m2 .  This  fact  provides  the  link  between  Gaussian  fields  and  Markov  processes  discussed 
in  the  Introduction. 

In  fact,  if  g°  is  the  Green  function  of  a  symmetric  Markov  process  A',  then  it  is  trivial 
that  g9  is  positive  definite,  so  there  exists  a  Slg8)  or  Ml(g9)  indexed  Gaussian  field  with 
covariance  kernel  g9 .  These  fields  are  called  the  fields  associated  with  A\  Although  not 
every  covariance  kernel  is  also  a  symmetric  Green  function,  this  is  so  often  the  case  that 
it  is  worthwhile  to  study  Gaussian  fields  from  the  viewpoint  of  Markov  processes. 

3b.  The  spectrum  and  Wiener-Ito  integrals. 

For  the  purposes  of  this  sub-section  we  shall  restrict  ourselves  to  the  parameter  space 
5,j.  We  call  a  random  field  stationary  on  Sd  if  $(/)=$(r2/)  for  all  /  €  S,i  and  /  £  RJ. 
where  Tzf(y)  =  f(y  +  x).  In  this  case  the  covariance  kernel  g  of  $  must  have  the  form 
g(x.y)  =  g{r  -  y).  and  so  we  write  it  as  a  function  on  RJ  only.  Furthermore,  (see.  e.g. 
Major  (1981)  for  details)  g  has  inverse  Fourier  transform  G,  which  is  a  <r-finite  measure  on 
R'1  satisfying  G\A)  =  G’(-A)  for  all  A  €  Bd  and 

(3.6)  J (l  +  ||A||)-QC(dA)  <  x 

for  some  ct  >  0,  (i.e.  G  is  tempered).  We  also  have 

£{*(/)*(/»)}=  [  /( A)h*(A)G’(dA).  f,h€Sd. 

J  R* 

where  ~  denotes  Fourier  transformation  and  *  complex  conjugation.  The  measure  G  is 
called  the  spectral  measure  of  $.  and  there  corresponds  to  it  a  random,  complex  valued. 
spectral  process  Z  =  Z(;.  defined  on  the  same  space  as  $,  such  that 

*(/)=  [  f(»ZG{d A). 


(3.7) 


f  €  Sd. 


The  importance  of  this  representation  for  us  is  that  not  only  does  it  give  a  simple  rep¬ 
resentation  of  stationary  Gaussian  fields  on  Sj.  but  it  also  leads  to  a  neat  characterisation 
of  all  their  C2  functionals. 

To  develop  this,  we  define  the  Fock  space  r<$  (or  TG)  of  the  field  $  as  the  Hilbert  space 

of  "Fock  columns"  f  =  (/0.  /i-  . . .).  where  /0  is  a  complex  number,  /„  =  /„  (Ai, _ An ). 

n  >  1,  complex  functions  of  At _ _  An  £  RJ,  satisfying 

(3-8a)  /*(A, . A»)  =  /*(-A, . -An), 

and  such  that  the  norm  | jfj)  admits  the  expansion 

(3.86)  ||f||2  =  l/o|2  +  f  >!)"*  /  •  •  •  /  |/*(A„  •  •  • ,  An)|2G(<fA1 ) . . .  G(dAn  )  <  oc. 

n  =  l  J  J 

It  is  well  known  (e.g.  Major  (1981))  that  there  exists  a  natural  unitary  mapping 
(of  Hilbert  spaces)  from  r<j>  to  C2[P)  defineable  via  multiple  Wiener-Ito  integrals,  taking 
f  £  T*  to  vPf  via 

(3.9)  :=/0  +  f>!p1  /  •  /  . An)ZG(dA1)...ZG(dAn). 

n  =  1  J  J 

with  H^fll2  =  ||f||2.  The  integrals  here  are  with  respect  to  the  spectral  measure  of  (3.7). 
and  do  not  charge  the  diagonal  At  =  A2  =  . . .  =  An. 

This  mapping  distinguishes  itself  from  all  other  unitary  mappings  T*  — *  C2{P)  by 
the  fact  that  the  Fock  columns  with  fm  —  0  for  m  >  n  correspond  to  polynomials 

p(Q{t j() _ ,Q(hm)),ht  £  Sd,  of  degree  n  in  the  values  of  the  field.  Other  important 

properties  of  this  representation  will  become  clearer  later. 

3c.  Wick  powers. 

One  way  to  get  away  from  a  purely  Gaussian  theory  of  random  processes  while  at 
the  same  time  retaining  as  much  as  possible  of  the  inherent  analytic  simplicity  of  this 
case  is  to  work  with  simple  functions,  such  as  powers,  of  them.  In  dealing  with  fields,  the 
corresponding  notion  is  that  of  a  Wick  power.  For  example,  if  $  is  defined  on  Sj  or  5(g) 


as  $(/)  =  /  f{x)$(x)dz.  then  we  would  like  to  define  a  k- th  power  :  :  that  would  art 

like  a  pointwise  k- th  power  of  i.e. 

(3.10)  :**(/)  •=  j  f(r)*k(x)dx. 

Of  couror,  if  (3.3)  is  ill-defined  then  so,  a  fortiori ,  is  the  above.  Nevertheless,  there  are  a 
number  of  ways  of  making  mathematics  out  of  (3.10). 

The  most  straightforward  would  be  to  consider 

(3.11)  lim  /  f(x)ik{d”)dx 

ft— oc  J 

as  a  candidate  for  :  $*(/)  where  . is  a  sequence  of  functions  in  S,i  or  S (3) 

converging  in  some  sense  to  the  delta  function  at  z.  Note  that  this  approach  also  works 
for  measure  indexed  fields.  Details  of  this  approach  are  implicit  in  Dynkin  (1984). 

As  opposed  to  the  ill-defined  (3.10)  and  the  somewhat  awkward  (3.11),  a  very  neat 
definition  of  Wick  powers  comes  by  identifying  &!:$*(/):  with  the  coefficient  of  tk  in  the 
formal  power  series  expansion  of  exp(/$(/)  -  £,,(/, /)].  That  is 

(3.12) 


_  V  L  **(/)  • 

k\ 

k=0 


Note  that  we  can  replace  /  here  with  /1  to  get  measure  indexed  Wick  powers. 

Although  (3.12)  is  neat  (and  useful)  it  sheds  little  light  ou  what  a  Wick  power  really 

is. 

Since  Wick  powers  belong  to  Z2{P)  iuey  also  have,  in  the  stationary  case,  a  Fork 
space  representation  in  terms  of  multiple  Wiener-Ito  integrals.  In  fact 

(3.13)  :**(/)•=  J  -  J  f(^+-+^)ZG(dXl)...Z(;(dXk). 

This  representation  does  not  carry  over  easily  to  either  the  non-stationary  or  the  measure- 
indexed  case.  Furthermore,  although  it  is  tidy,  in  no  way  is  it  obvious  why  a  spectral 
integral  of  this  form  should  correspond  to  the  pointwise  power  (3.10). 

It  turns  out  that  Wick  powers  have  a  veiy  natural  representation  in  terms  of  the 
Markov  path  intersections  mentioned  in  the  Introduction,  and  this  sheds  light  both  on 
their  structure  and  their  uses,  some  of  which  we  shall  see  in  Section  5  below. 
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3d.  Additive  functionals. 

The  last  concept  that  we  shall  need  for  the  moment  relating’  to  Gaussian  fields  is  that 
of  additive  functionals.  To  define  these,  let  At-,  denote  set  of  all  closed  intervals  in  R1 
of  the  form  [a.b]  =  {x  :  a,  <  x,  <  b,}.  (a  <  b).  Let  A  denote  the  ring  of  finite  unions 
Ai  U  . . .  U  An  of  elements  of  A0  with  pairw  ise  disjoint  interiors,  i.e.  A,  D  A}  =  0.  «  ^  j. 

Write  B*  for  the  <r-algebra  generated  by  i.e.  B*  =  <r {♦(/)}  if  the  field  is  Sj  or  S|<;) 
indexed,  or  B#  =  <r{$(ji)}  if  the  field  is  indexed.  Also,  for  open  0  C  RJ  define  B<t>,u 

to  be  the  <r-algebra  generated  by  $(/)  with  supp /  C  0,  or  $(**)  with  suppji  C  0.  as  the 
case  may  be.  Then  if  to  every  A  €  A  there  corresponds  a  random  variable  €  £2(P). 
measurable  with  respect  to  B<$,  we  call  ♦  :  A  —  £2(P)  a  tquare  integrable  functional  of  $. 

Now,  for  closed  A  C  Rd,  set  B*,^  =  floo  &*.ai  ,  where  A ‘  is  the  t-neighbourhood  of 
A.  Then  there  are  three  properties  of  square  integrable  functionals  that  will  be  of  interest 
to  us: 

I.  Locality.  For  each  A  £  A  the  random  variable  VtA  is  measurable  with  respect  to 
B<p,  a  ■ 

II.  Additivity.  *AUB  =  9 A  +  P- a.s  for  ail  A,B  €  A  with  °A  n  ‘B  =  0. 

III.  Stationarlty.  =  U^9A,  where  Ux,x  E  Rd,  are  the  unitary  translation 

operators  on  £2(P)  corresponding  to  the  translations  Ujf(y)  =  f[y-  x)  if  /  €  Sd  or  S(ff), 
or  to  the  translations  UX^(A)  =  n[A  -  x)  if  /i  €  M1^). 

The  central  interest  in  these  three  properties  lies  in  the  fact  that  local,  additive, 
stationary  functionals  can  be  applied,  in  a  relatively  simple  fashion,  to  stationary,  Markov, 
Gaussian  fields  to  yield  stationary,  Markov  non-Gaussian  fields.  The  prescription  is  merely 
to  change  from  the  Gaussian  measure  P  to  a  new  measure  exp{¥R<}  •  P.  For  details  see, 
for  example,  Rozanov  (1982). 

To  study  these  three  properties  in  the  case  of  stationary  Gaussian  fields,  Dobrushin 
and  Kel'bert  ( 1983a, b),  noted  that  since  each  ♦a  E  £2(P),  it  must  have  a  Wiener-Ito 
integral  representation  of  the  form  (3.9).  If  we  assume,  for  no  loss  of  generality,  that 
£•{*.,  }  e  0.  then  we  call  the  Fock  column  - )  in  this  representation  the  spectral 


representation  of  Dobrushin  and  Kel'hert  showed  that  it  is  possible  to  read  off  the 
above  properties  from  the  spectral  representation  of  We  shall  investigate  their  result  in 
Section  5,  via  completely  different  methods,  and  give  it  an  intuitive  explanation  not  linked 
to  the  stationary  situation. 

We  now  have  (somewhat  more  than)  enough  background  to  look  at  our  first  new 


result. 


4.  The  Limit  Theorem  and  its  Consequences. 

We  now  return  to  the  setting  of  the  Introduction,  with  one  small  change.  Let  (ft.  7 .  P) 
be  a  probability  space  rich  enough  to  support  the  following: 

(a)  a  Poisson  random  variable  .V  =  A’,\  with  mean  A  >  0. 

(b)  an  infinite  sequence  of  independent  right  Markov  processes  A’i  (/),  X2(/) . t  >  0. 

taking  values  in  and  with  common  stationary*,  symmetric  transition  density.  pt{ x,  y ). 

(c)  an  infinite  sequence  of  independent  variables  oq.tr 2 _ _  with  P(cr,  =  +1)  = 

P(<r,  =  -1)  =  i. 

The  two  infinite  sequences  and  A'  are  all  independent  of  one  another. 

Fix  t)  >  0.  and  let  ge{ /,  y)  be  the  Green  function  of  the  A, .  Fix  A  >  0.  take  f  €  S{ge ), 
and  define 

,  NU) 

(4.1)  **Klf):=[9/\)2 

i  =  L 

This  is  the  type  of  sum  we  looked  at  in  (1.5)  as  being  equivalent,  in  that  case,  to  the 
ill-defined  sum  of  integrals  of  local  time.  Note  that,  by  varying  /  in  (4.1).  becomes  a 
random  field  on  Sobolev  space.  This  field  has.  as  A  — *  oc.  a  Gaussian  limit. 

Theorem  4.1.  As  A  —  :c  the  random  held  converges,  in  the  sense  of  convergence  of 
finite  dimensional  distributions ,  to  the  zero  mean  Gaussian  Geld  on  S((7#)  with  covariance 
kernel  g9 . 

This  result  is  actually  a  special  case  of  the  more  general  Theorem  4.2  below,  but  it  is 
worthwhile  to  state  it  separately  because  of  its  simplicity.  (Both  will  be  proven  in  Sections 
7  and  8.)  The  result  is  also  true  if,  throughout,  we  replace  the  Sobolev  space  S(05)  by  the 
Schwartz  space  $d,  and  impose  the  extra  condition,  on  g8 ,  that 

f[jc)g0[x,  y)f(y)didy  <  oc.  for  all  /  G  Sd- 

Th*'  extension  to  fields  defined  on  the  space  of  measures  M1  (ge)  is  not  so  easy  to  formulate. 
Muce  we  need  to  work  a  little  harder  defining  what  we  mean  by  $*x(//)  (The  Gaussian 


limit  process  is  easy  to  identify,  however.)  To  treat  the  measure  case  we  now  move  to  a 
far  more  general  formulation,  although,  for  reasons  that  we  shall  note  later,  we  do  not  yet 
want  to  treat  the  most  general  formulation.  Both  allow  us  to  consider  not  just  CLT's  for 
the  Markov  A’,,  themselves  but  also  for  the  interaction  of  k  such  processes. 

Fix  k  >  1,  6  >  0,  and  consider  the  set  of  measures  7  =  ~f(dri . dik )  on  that 

belong  to  M9k  =  Mk{ge)  and  satisfy  the  symmetry'  condition 

(4.2)  7(-4, 4*)  =  7(.4,j . A,J 

for  each  permutation  (i| . fy)  of  (1 . k ).  and  A,  €  Bd .  We  denote  this  collection  by 

Kymm  or  Msymmid8)- 

By  Theorem  2.1.  there  exists  an  additive  functional  F 9  corresponding  to  each  7  6 
defined  on  k  of  the  X,  at  a  time.  Consider  firstly  the  case  k  =  1,  and  define 

,  M-M 

(4.3)  <(7):=(0/A)2  £>,ft(A\). 

>=1 

Then  $ex  can  be  considered  to  be  a  random  field  on  and  an  analogue  of  Theorem 

4.1  holds  for  it. 

Theorem  4.2.  As  A  —  oc  the  random  field  $ex  of  (4.3)  converges,  in  the  sense  of  conver¬ 
gence  of  finite  dimensional  distributions,  to  the  zero  mean  Gaussian  held  on  with 

covariance  kernel  ge . 

Thus,  between  Theorems  4.1  and  4.2  we  now  know  how  to  approximate  both  function 
and  measure  indexed  fields  via  Markov  paths.  (Note  that  for  a  y  6  with  density 

/  €  5 (/),  we  have  •»  s  K  (/),  and  both  theorems  identify  the  same  limit  random 
variable.) 

We  would  now  like  to  approximate  elements  in  the  Z2  space  of  a  Gaussian  field.  Since 
our  discussion  in  the  previous  section  relied  on  the  Fock  space  representation  of  the  Z 2 
space,  and  this  required  stationarity  of  the  field,  we  shall  assume  for  the  moment  that 
Pt(jr.y)  =  pt[r  -  y).  or.  equivalently,  0H(Ay)  =  g?(r  -  y).  This  implies  (Blumenthal  and 


Getoor  (1968),  p.  17)  that  the  A',  have  independent  increments.  We  shall  lift  this  condition 
later,  as  we  shall  also  do  with  the  following  conditions: 

(a)  /  /  /M/(*  ~  y)f{y)dxdy  <  co  for  all  /  6  $dt 

(b)  g6  =  C»  ,  for  some  <r-finite,  symmetric,  tempered  measure  Ge  on  Rd, 

(c)  e  has  a  Fourier  transform  7  which  exists  as  a  regular  (i.e.  non-generalised) 

complex  valued  function. 

We  write  the  set  of  7  satisfying  (c)  as 

It  follows  from  (a)-(c),  and  a  generalised  Parseval-type  equality  (c.f.  the  generalised 
Poisson  summation  formula  of  Argbright  and  Gil  de  Lamadrid  (1972))  that  7  6  where 

:=  {/  =  /( Ai . A *),  A,  6  Rd;  ||/|&x...xG  <  00  and 

/(A, . Afc)  =  /*(-A, . -Afc)  =  /(A, A, J) 

where 

ll/llox...xo  =  /•••/ l/(Ai . A*  )|aG(dA, ) . . .  G(rfAfc ). 

Consequently,  if  we  write  Qe  for  the  stationary  zero  mean  Gaussian  field  on  Sj  with  spectral 
distribution  function  G9 ,  then  the  functional 

(4.4)  *’:=F /•  ••  /l(A,,...,A*)2c.(<iA,)...2G-(<fAt) 

is,  according  to  Section  3b,  well  defined  as  an  element  of  the  £}  space  of  Q9.  Furthermore, 
there  corresponds  to  7  a  functional  F9  of  k  independent  Xt.  The  link  between  7,  ^  and 
F9  is  via  the  sum 

ro  ,  N(X)  <  k, 

(4.5)  *1(7):=  W0/A)*/2  E---E  ,N(X)>k , 

for  which  we  have 

Theorem  4.3.  Under  tbe  above  conditions,  the  pair  ($*(/),  ^(7))  converges  in  distri¬ 
bution.  as  A  —  00,  to  the  pair  {**(/),  (7))- 

Remark.  Note  that,  although  we  have  not  stated  Theorem  4.3  in  such  terms,  by  varying  / 
over  Sd  and  7  over  *M.9Symm-  we  could  write  the  theorem  as  being  about  convergence  of  the 


vector  processes  {$^(-).  ¥*(-))  on  5,;  x  *Ms'vfcmm  to  the  vector  process  ($*(•).  4,!,(-|).  with, 
as  before,  convergence  being  in  terms  of  finite  dimensional  distributions.  This  extended 
form  of  the  theorem  is  an  immediate  consequence  of  the  linearity  of  all  the  fields  on  their 
respective  parameter  spaces. 

Thus,  modulo  the  present  technical  restrictions  on  ge  and  7.  we  have  found  a  way  to 
represent  Gaussian  fields  via  Markov  processes,  and  functionals  of  the  former  by  functionals 
of  the  latter.  The  technical  conditions  will  all  be  lifted  in  Section  8,  but  their  imposition 
at  this  stage  makes  the  didactic  discussions  of  the  following  section  somewhat  easier  to 
follow.  These  discussions  hinge  on  using  Theorem  4.3  to  think  of  Gaussian  fields  and  their 
functionals  in  terms  of  approximating  sums  of  functionals  of  Markov  processes.  However, 
because  of  the  restrictions  we  have  placed  on  y,  we  cannot  consider  all  functionals  of  either 
type.  Nevertheless,  the  next  two  lemmas  show  that  we  can,  via  these  restricted  classes  of 
functionals,  get  arbitrarily  close  to  any  functional  of  interest  (in  the  stationary  case).  The 
lemmas  and  the  discussion  following  them  can  be  skipped  on  the  first  reading. 

Lemma  4.1.  *Mf>sktnm  is  everywhere  dense  in  Mgkmm. 

Lemma  4.2.  The  set  of  linear  combinations  of  f  with  inverse  Fourier  transform  in 

X5yram (g0)  is  everywhere  dense  in 

The  consequences  of  these  lemmas,  which  we  shall  establish  in  a  moment,  are  immedi¬ 
ate,  for,  together  with  Theorem  4.3,  they  tell  us  that  any  Gaussian  field  on  S(ge),  together 
with  any  of  its  functionals,  has  an  approximation  via  Markov  processes,  (assuming,  of 
course,  that  the  covariance  kernel  is  also  the  Green  function  of  a  right  Markov  process). 
Proof  of  Lemma  4.1.  Consider  7  G  Mesykmm.  We  shall  build  a  sequence  -^.6  = 

(^i . ).  with  G  *Mes'ykmm  such  that  7^  —  7  as  j£|  —  0.  Convergence  is  in  terms  of 

the  metric  defined  by  (2.10),  i.e.  via  the  inner  product  (•,  ■)#. 

Firstly  note  that  by  a  standard  truncation  argument  the  finite  7  G  are  dense 

in  Mea'kmm,  and  so  we  may  assume  7  to  be  finite.  Now  consider  the  absolutely  continuous 


approximation  to  7  used  in  constructing  additive  functionals  in  §2b  and  given  by 

7 g(dx)  ■■=  be^g(x)dx 

=  j  e~elslpB(x.y)^(dy)dx. 


But  h*  ,  €  £‘(R'ik)-  since 

i.O 


J  |6®tf(x)|dx  =  J  dx  J  i[dy)e  8|tf|p^(x.y; 

=  e-9^1  J  7(rfy)  J  P$(x.y)dx 


<  oo. 


Thus  6**  «  has  a  Fourier  transform  which  is  a  continuous,  non-generalised,  function  and  so 

7.0 

70  €  Furthermore,  it  is  trivial  to  check  that  (~jg  -  7,  7$  -  "))e  —  0  as  |6j  — ►  0,  so 

we  are  done. 

Proof  of  Lemma  4.2. 

Let  5j  denote  complex  Schwartz,  space,  so  that  f  €  Sd  is  of  the  form  /  =  / 1  +  1/2. 
fi.h  &  Sd ■  Assume,  for  the  moment,  that  k  =  1,  and  note  that  C  £2(G<?).  Then 

since  Sd  is  dense  in  £2(Gd)  it  is  easy  to  see  that 


cc 

"dysymm 


:={/eSf:n*)  =  /(-*)> 


is  also  dense  in  .  Thus,  for  the  case  k  =  1,  it  will  suffice  to  show  that  an  f  €  Sd iSJ/mm 
can  be  written  as  the  sum  of  the  Fourier  transforms  of  two  measures  in  At£ymm  (A 

To  see  this,  note  that  the  mapping  f  —*  J  is  an  invertible,  bicontinuous  transformation 
from  Scd  to  itself.  Thus  if  /  €  Scd ,symm,  then  F~lf=f€S cd  where  F~l  represents  inverse 
Fourier  transform.  Furthermore,  /  is  real,  since 


M  ^  /  '‘*/(*H* 


=  (/'PI  =  /(-■.)! 


Writing  h+  —  raa x(h,  0).  h  =  -  min(/o  0).  we  have,  for  /  e  Sd  symm. 

(ir-,/)(A)  =  (/(A))+-(/(A))' 

with  both  /+  and  /~  real.  To  obtain  measures,  write  ■)  +  (dA)  =  /+(A)dA.  ~J{d\)  = 
f~{\)d\.  Note  that  by  the  generalised  Poisson  summation  formula, 

oo  >  j \f\2G{dX) 

=  J  fi*)g9{*  -  y)f{y)dxdy 

=  J  f+(x)ge(x  -  y)f+(y)dxdy  +  J  f~  (x)ge  (/ -  y)f~  (y)dxdy 
-2  J  f+(x)g9(x  -  y)f~(y)dxdy. 

Since  /  6  5<*,  and  so  is  arbitrarily  small  outside  large  compacts,  each  of  the  last  three 
integrals  must  be  finite.  Consequently  both  7  +  and  ~ij  have  finite  (•,  • )8  norm  since, 
for  example,  {'/,  =  /  f+{x)g6{x  -  y)f+{y)dxdy.  Thus  both  7+  and  7J  are  in 

Kym  m(g° )  and  we  are  done. 

The  proof  for  general  k  >  1  follows  by  noting  that  the  elements  of  x  . . .  x 
are  dense  in  and  that  for  /  €  of  the  form  /(*!,...,**)  =  /1  ( x i ) . . .  /*  ( Xk )  with 
/,(/,)  6  "Hq,  we  have  /  =  A  ...  A.  The  result  for  k  —  1  can  then  be  applied  to  write  /  as 
a  linear  combination  of  measures  in  .M* 


5.  Applications,  I. 

In  this,  aud  the  following  section,  we  return  to  what  originally  motivated  us  -  obtain¬ 
ing  insight  into  random  fields  via  approximating  sums  of  Markov  processes.  The  results 
range  from  intuitive,  non-rigorous  arguments  as  in  our  discussion  of  the  Markov  property 
of  random  fields,  (which  generate  no  new  results)  to  rigorous  arguments,  such  as  the  dis¬ 
cussion  of  renormalisation  in  the  following  section  where  the  results  are  essentially  new. 
We  start  with  what  is  basically  an  example  for  Theorem  4.3.  the  construction  of  Wick 
powers. 

5a.  Wick  powers. 

In  )}3c  we  gave  a  sequence  of  equivalent  definitions  of  Wick  powers  (see  also  §8a), 
including  (3.13)  which,  for  /  €  Sj.  set 

(5.1)  :**(/)  :=  J  ...  J  f\\l+...  +  \k)ZG,(d\l)...ZG,(d\k). 

for  $  with  covariance  kernel  ge .  Consider  the  measure 

(5.2)  ';(<*/ 1 . dxk)  ■-  /Ui^Ui  -  *i)-- ~  *k)dii  ...dik, 

where  d  is  the  Dirac  delta  function.  (Formally  7  is  merely  the  measure  concentrated  on 
the  diagonal  xx  =  . . .  =  xk  in  (lRd)fc,  with  density  /Ui).)  The  Fourier  transform  of  7  is 
given  by 

7  ( A 1 , . . . ,  Afc)  =  J  e'~x>x,f{xf)d{  xi  -  x2)...b(xk^1  -  xk)dxx  ...dz  k 

=  j  e'xZX’f(x)dx 
=  /( Ai  +  . . .  +  A* ). 

Thus,  comparing  with  (5.1).  we  have 

:»*(/):=  j  ■■  j  ;(A, . \k)ZGi(d\i ) . .  .  ZG*(d\k), 

which  is  precisely  k  1  times  the  ^  of  (4.4)  and  Theorem  4.3. 


Consequently,  recalling  the  properties  of  measures  like  7  discussed  in  §2c.  we  have  a 
representation  of  Wick  powers  based  on  intersections  of  paths  of  Markov  particles  the 
A;-th  power  coming  from  intersections  of  k  particle  paths. 

This  representation  will  work  as  long  as  7  6  Clearly,  it  is  only  the  integra- 

bility  of  7  that  is  of  importance  here.  If  our  particles  execute  Brownian  motion,  so  that 
the  limit  Gaussian  process  is  the  free  field  of  §3a.  it  is  easy  to  check  that 


(5.3) 

d  =  2  => 

•v  g  *Ma-k 
•  * %symm 

for  all  k. 

(5.4) 

d  =  3  =» 

ci  *k(Q'k 

1  c  J  xgymm 

only  for  k  = 

(5.5) 

d  >  4  => 

1  J  xsy mm 

only  for  k  = 

That  is,  this  technique  would  seem  to  work  only  in  very  restricted  cases.  These  are 
precisely  the  cases  in  which  k  independent  Brownian  motions  in  Rd  have  intersecting  paths 
(c.f.  Dvoretsky  et.  al.  (1954).) 

Nevertheless,  if  we  recall  that  Gaussian  fields  also  go.  in  Physics,  under  the  name  of 
Euclidean  quantum  fields,  then  we  can  verify  (e.g.  Glimm  and  Jaffe  (1982))  that  Euclidean 
field  theory  turns  out  to  be  interesting  only  in  the  cases  in  which  our  limit  theorem  applies 
to  Wick  powers,  so  we  have  not  done  too  badly. 

Before  leaving  this  example,  we  note  that  it  is  essentially  this  problem  with  d  =  2  and 
A'  =  Brownian  motion  that  the  two  papers  of  Wolpert  (1978a,  b)  are  dedicated  to.  For 
more  information  on  the  functionals  related  to  F *  with  7  given  by  (5.2)  with  /  =  1.4.  see 
the  papers  on  intersection  local  time  by  Geman,  Horowitz  and  Rosen  (1984)  and  Rosen 
(1985).  For  more  serious  applications  of  this  to  Physics,  see  Aizenmann  (1982).  The 
original  idea  of  using  intersection  local  time  to  study  Wrick  powers  comes  from  the  seminal 
paper  of  Symanzik  (1969). 

5b.  Functionals  of  Gaussian  fields. 

In  §3d  we  introduced  the  notion  of  functionals  of  Gaussian  fields,  along  with  three  pos¬ 
sible  properties:  additivity,  locality,  and  stationarity.  Dobrusbin  and  Kel'bert  ( 1983a. b) 
studied  these  from  the  point  of  view  of  Fock  space,  which  they  used  to  write  a  general 
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functional  via  its  spectral  representation.  Consider  the  k- th  term  in  such  a  representation, 
defined  by 

(5.G)  -  —  j  ...  j  vk\  Ai , . . . ,  A*  )Zc(rfA| ) . . .  Za(d\k  )• 

To  study  the  properties  of  tyk  as  A  varies  over  unions  of  rectangles.  A.  let  irk  be  the 
natural  embedding  of  the  space  W*.  into  the  space  S'ik  of  generalised  functions  and  define 
the  mapping  Jfk  :  A  —  S'lk  by 

=  **U-£). 

Then  Dobrushin  and  Kel'bert  have  shown  the  following: 

I:  'Pfc  is  additive  iff  nk  is  additive,  in  the  sense  that 

(5.7)  ?FfcU,  U  A2)  =  5Ffc(A,)-f  jffc(A2). 

for  disjoint  .4i,.42  6  A. 

II;  is  local  and  additive  iff  Jfk  is  additive  and  diagonal,  in  the  sense  that 

(5.S)  supp  F-1  |f*,(.4))  C  diagAfc  :=  {(/i . xk)  :  =  ...  =  /*  e  .4}, 

where  F~x.  as  before,  represents  inverse  Fourier  transform. 

Ill;  i'fc  is  stationary  iff  F~lWk  is  invariant  in  the  sense  that 

(5.9)  Oa+A/) 

where  0A  6  S^k  is  given  by  F_1F*,-(A),  /  €  Sdk ■  *  €  diagIR1*)*.  and  Uz  is  the  shift 
operator  of  §3d. 

Note  two  things.  Firstly,  since  everything  is  based  on  the  Fock  space  representation, 
the  above  results  are  valid  only  when  the  underlying  (iaussian  field  is  itself  stationary. 
Secondly,  although  these  results  are  precise,  they  have  little  intuitive  appeal.  In  particular, 
we  ask.  why  should  the  inverse  Fourier  transform  of  the  kernel  v*  figure  so  prominently 
in  the  conditions  ’  The  answer  to  this  is  clear  once  we  approach  the  problem  from  the 
viewpoint  of  approximating  fields  and  functionals  by  Markov  paths. 


We  now  need  to  assume  that  the  underlying  field  has  covariance  kernel  g* .  which,  as 
usual  is  a  Green  function.  Consider  the  sum  ¥^(7)  of  (4.5)  with  7  =  7*.  where  7^  =  l'  *  . 
(Note,  to  ensure  that  v*  has  an  inverse  Fourier  transform  in  m ,  as  7^  must  he.  we 

really  need  to  assume  this.  However,  since  Lemma  4.2  implies  that  such  v  are  dense  in  M(k, „ , 
we  lose  no  real  loss  of  generality  by  tacitly  doing  so.)  By  Theorem  4.3,  ¥^(.4)  :=  ^(7^) 
converges  in  distribution  to  the  of  (5.6).  and,  indeed,  does  so  as  a  process  over  .4  6  A. 

The  first  thing  we  note  is  that  the  measures  7^  now  become  the  natural  parameters 
of  the  problem,  and  so  one  can  see  where  the  inverse  Fourier  transforms  of  (5.8)  and 

(5.9)  might  come  from.  In  fact,  even  (5.7)  could  be  written  in  terms  of  inverse  Fourier 
transforms,  since  it  is  clearly  equivalent  to 

(5.10)  F~lnk(Al  U  A2)  =  F_1  jf(Ai )  +  F-1jrfc(A2). 

Secondly,  since  F-1iFfc(A)  =  7^,  and  the  7^  and  F^A  are  easy  to  understand,  some 
intuition  may  now  be  forthcoming.  We  develop  this  intuition  by  studying  additivity,  local¬ 
ity.  and  stationarity  not  for  the  functional  of  the  Gaussian  field,  but  for  the  functional 
V\(7i4)  of  the  Markov  processes  which  ultimately  converges  to  ♦£.  Note  that  we  cannot 
rederive  the  results  of  Dobrushin  and  Kel’bert  this  way,  but  can  only  provide  a  way  of 
understanding  them. 

a.  Additivity.  It  is  immediate  from  the  definition  of  ^(7^ )  and  F®„  that  the  former 
is  additive.  In  fact,  our  construction  has  always  been  such  that  we  can  build  only  additive 
functionals,  since  7^  is  required,  a  priori ,  to  be  a  measure.  Condition  (5.7)  (or  (5.10)) 
simply  allows  our  construction  to  be  started. 

b.  Locality.  Here  we  shall  have  to  define  a  slightly  amended  definition  of  locality  more 
suited  to  the  Markov  path  setting.  For  any  open  O  C  R*.  we  let  So  denote  the  <r-algebra 
generated  by  the  .Y,  while  they  take  values  in  O.  This  notion,  which  we  shall  call  the 
excursion  cr-field,  can  be  made  precise.  (See.  for  example  Jacobs  (1978).  Kaspi  (1985)). 
We  now  define  locality  as  in  §3d,  but  with  respect  to  these  (T-algebras. 

Assume  (5.8),  i.e.  supp  7^  C  diag.4*.  It  is  then  clear  from  the  definition  of  the  Fri, 
that  they,  and  so  ^(7^).  are  local  in  the  new  sense.  Furthermore,  it  is  easy  to  see  that 


this  follows  from  supp  7^  C  Ak.  the  additional  restriction  to  the  diagonal  being  seemingly 
superfluous.  This  restriction  arises  in  a  natural  fashion  from  the  additivity  of  as  in 
Dobrushin  and  Kel  bert  (1983a).  Using,  then,  what  we  know  about  functionals  of  Markov 
processes,  we  see  that  the  diagonality  condition  implies  that  the  k- th  order  functional  of 
the  Gaussian  field  is  approximated  by  a  functional  built  on  the  intersections  of  k  particles 
at  a  time.  This  is  the  full  import  of  (5.8). 

Thus,  with  sufficient  motivation,  we  could  build  an  approximation  to  the  whole  of 
Fock  space,  the  k- th  part  of  this  based  on  functionals  of  k  particles.  The  reader  familiar 
with  quantum  field  theory  will  see  an  extremely  pleasant  semantic  coincidence  here,  for 
the  quantum  theory  of  k  interacting  particles  is  also  the  theory  of  the  fc-th  part  of  Fock 
space.  Needless  to  say,  our  particles  and  those  of  quantum  field  theory  have  no  physical 
link  semantically  however,  the  coincidence  is  intriguing. 

c.  Stationarity.  From  their  construction  the  functionals  ^(7^)  are  stationary  in  the 


sense  that 


Ux<bfc4)  =  <(U,7fc4). 


where  Vx  is  the  unitary  translation  operator  on  the  £2  space  of  the  Markov  processes 
corresponding  to  shifting  all  starting  points  of  the  paths  by  -r.  and  Uz  is  defined  by 


/  Ul‘ 
J  B 


(dr, . dxk 


= L 1 


(dr, . drfc). 


Consequently,  if  we  want  them  to  be  stationary  in  the  sense  of  Dobrushin  and  Kel  bert 


r>r  =  *!b*4+x). 

(c.f.  §3d )  then  we  require  that  74+x  =  f ',7*  ■  But  this  is  precisely  condition  (5.9).  Thus 
we  see  that  this  condition,  too.  arises  as  a  condition  ou  the  paths  of  our  Markov  particles. 

5c.  The  Markov  property  for  Gaussian  fields. 

One  of  the  most  interesting,  and  most  deeply  studied  properties  of  Gaussian  field*  is 
that  of  their  Markovianess.  Unlike  the  case  for  regular  processes  m  univariate  time,  the 


Markov  property  for  measure  of  function  indexed  fields  is  not  particularly  useful.  Nev¬ 
ertheless.  it  is  probabilistically  intriguing.  Furthermore,  it  is  of  considerable  importance 
in  Physics,  since  in  the  transformation  of  regular  field  theory  to  Euclidean  field  theory 
(i.e.  to  Gaussian  fields)  one  of  the  axioms  of  the  former  transforms  to  an  assumption  of 
Markovianess.  (c.f.  Nelson  (1973)).  Consequently,  it  is  only  the  Markov  Gaussian  fields 
that  are  of  interest  to  physicists. 

There  are  many  different  definitions  of  Markovianess.  and  they  are  not  all  equivalent. 
Albeverio  and  Hoegh-Krohn  (1984)  give  a  detailed  discussion,  but  we  shall  follow  the  treat¬ 
ment  of  Dynkin  (1980).  Firstly,  we  restrict  our  attention  to  measure-indexed  (Xs -1)  fields. 
We  say  that  a  field  $  is  Markov  with  respect  to  two  sets  A,  B  C  R*  if  the  tr-algebras  /<*>,. * 
and  7*.b  are  conditionally  independent  given  7*. ,4ns-  where  J^  A  :=  <r{$(/i)  :  supp  n  C 
.4).  (Note  that  7<p_A  is  not  the  same  as  the  £<*>  4  of  §3d.)  One  easily  extends  this  to  a 
general  Markov  property  by  demanding  that  for  all  B  C  Rd.  with  smooth  enough  bound¬ 
ary  dB  and  complement  Bc .  <&  is  Markov  with  respect  to  B  and  Bc .  It  is  clear  therefore 
why  it  is  more  natural  to  work  with  measure  indexed  fields  would  be  empty  in 

the  function  indexed  case,  for  we  could  generally  find  no  /  G  Su  with  support  in  OB.  (A 
way  around  t hi-,  is  to  base  the  Markov  property  on  £$,  4  instead  of  7<*,.. 4.  This  leads  to 
a  different  theory  which,  when  specialised  to  simple  processes  on  R1.  for  example,  is  not 
the  usual  Markov  theory,  c.f.  Rozanov  (1982).  Wong  and  Zakai  (1986).) 

If  we  limit  our  attention,  as  usual,  to  Gaussian  fields  with  covariance  kernel  the  Green 
function  of  a  Markov  process,  then  Dynkin  (1980)  has  a  very  elegant  characterisation  of 
Markovianess 

Characterif  atlon.  $  has  the  Markov  property  on  all  sets  A.B  for  which  it  is  impossible 
for  thi  associated  Markov  process  to  reach  A  from  B  without  crossing  A  H  B.  (This  result 
also  has  a  ric  r.-  precise  statement  see  both  Dynkin  (1980)  and  \tkmson  (1983).) 

V\ e  would  like  to  understand  this  result  from  our  viewpoint.  Iu  fact,  it  will  suffice  if 
we  i-au  see  why  the  condition  of  the  conjecture  implies  the  are  Markovian.  This  will 
not.  .if  course  imply  the  Markovianess  of  <l>  itself,  nor  can  it  ever  hope  to  (the  convergence 
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<>f  $ '  to  $  is  in  distribution  only,  and  th<‘  Markov  property  is  in  essence  an  a.s.  statement) 
but  it  does  provide  insight. 

Fix  A.  and.  for  the  sake  of  intuition,  fix  A’(A)  =  A  =  ri  in  the  sum  (4.3)  defining  $^(p). 
which  now  becomes 

n 

(5.12)  <(/i)  =  W«),/2 

t  —  1 

Now  we  argue  more  or  less  heuristically.  If  the  A',  have  regular  enough  sample  paths, 
then  we  claim  that  /($f(,.4)  is  identical  to  the  excursion  tr-algebra  generated  by  X  = 

(A’ i . A",,)  while  the  A',  are  in  .4.  Then,  from  rigorous  results  of  Jacobs  (1978)  and 

Kaspi  (1985)  one  can  see  that,  again,  requiring  regular  paths.  7($^..4)  and  J($®,AC) 
will  be  conditionally  independent  given  7($^.d.4).  But  to  apply  these  results  the  balance 
between  regularity  of  the  sample  paths  of  the  A,  and  dA  must  be  such  that,  at  very 

least,  one  can  define  an  additive  functional  of  (A'i . A'n)  living  on  dA.  Given  the  other 

assumptions  that  we  have  placed  on  the  A,,  this  means  that  we  require  the  A,  to  be 
continuous,  and  so  we  have  a  heuristic  version  of  Dynkin's  characterisation  result. 
Actually,  all  of  this  is  easiest  to  see  when  d  =  2.  n  =  1,  and  we  draw  a  picture. 


Figure  5.1.  Markov  paths. 

In  the  leftmost  case  example  of  Figure  5.1  the  Markov  path  is  continuous,  and  it  is 
clear  that,  given  all  the  hitting  points  of  the  path  on  3D.  each  excursion,  between  hitting 
points  is  independent  of  all  the  others.  Consequently,  all  the  excursions  contained  in  B 


ar<*  conditionally  independent  of  those  in  D  .  and  from  this  follows  the  Markov  property 
of  In  the  rightmost  case,  the  path  may  never  hit  dB.  and  so  there  is  no  information 
in  to  even  attempt  to  start  conditioning. 

Finally,  we  note  that  it  is  tempting  to  try  to  push  these  arguments  further,  to  obtain 
an  expression  for  £{$*(/i)  |  with  supp  n  6  Bc .  as  Dynkin  (1980)  has  done,  iu 

terms  of  the  first  hitting  time  of  .V,  on  dB ,  conditioning  on  starting  -V,  according  to  //.  It 
is  clear  how  the  excursion  processes  associated  with  A',,  along  with  time  reversibility,  need 
to  be  used  here,  but  the  calculations  involved  seem  so  complex  that  nothing  is  gained  via 
this  approach. 

5d.  The  free  field  again. 

It  is  instructive  to  close  this  section  with  an  almost  trivial  comment  on  the  free  field 
of  (3.5),  which  holds  a  uniquely  central  role  in  both  Euclidean  field  theory  and  Gaussian 
fields,  as  the  only  stationary  Markov  field.  Our  question  is  -  what,  from  our  point  of  view, 
makes  this  process  so  special?  We  know  that  to  build  a  Gaussian  field  we  need  Markov 
paths.  If  we  require  the  field  to  be  stationary  then,  as  we  have  already  noted,  the  Green 
function  must  satisfy  y)  =  g*[i  -  y).  This,  however,  forces  the  Markov  processes  to 
have  stationary,  independent  increments.  If  we  furthermore  demand  Markovianess  for  the 
field,  then  we  demand  continuity  for  the  paths  of  the  Markov  particles.  There  is  only  one 
continuous  process  with  stationary,  independent  increments  Brownian  motion.  Thus  we 
have  only  one  way  to  build  the  field  from  Markov  paths.  The  resulting  field,  as  we  have 
already  seen,  is  the  free  field. 


6.  Applications  II:  Renormalisation 


III  the  previous  section  we  used  the  limit  theorems  of  Section  4  primarily  at  a  heuristic 
level  to  explain  a  number  of  existing  results.  In  this  section  we  shall  use  these  theorems 
in  a  fully  rigourous  fashion  to  motivate,  and  then  prove,  new  results.  The  central  theme 
of  this  section  comes  from  the  idea  of  renormalisation  techniques,  which  we  shall  define 
formally  soon. 

6a.  Self-similarity  for  fields. 

In  Probability  Theory,  renormalisatiou  is  generally  (but  incorrectly)  considered  to  be 
synonymous  with  the  notion  of  self-similarity.  To  make  our  lives  notationally  easier,  we 
shall  restrict  our  attention  iu  this  section  to  random  fields  on  Sd,  in  which  case  we  call  a 
field  $  selJ-*imUa.T  of  order  n  iff  for  every  t/  >  0  and  /  G  $d 

(6.1)  $(/)="*(/):=*(*/) 

where 

(6.2)  ”/{*)  :=  t/-° /(•*•/»>)• 

This  is  clearly  a  scaling  phenomenon,  and  in  the  Gaussian  situation  Dobrushin  (1979)  has 
show  n 


Theorem  6.1.  Let  $  be  a  stationary ,  zero  mean,  Gaussian  random  field  with  spectral 
measure  G.  Then  $  is  self-similar  of  order  a  iff 

(6.3)  <7(.4)  -  rr'2o(r(nA),  .4  £  5(RJ\{0}),  t)  >  0. 


As  neat  as  this  result  is.  there  are  two  problems  with  it.  Firstly,  it  applies  only  to 
stationary  fields.  Since  both  the  statement  and  proof  rely  on  spectral  methods,  there  is 
no  clear  extension  to  the  non-stationary  case.  More  importantly,  however,  is  that  it  deals 
with  far  too  narrow  a  notion  of  self-similarity.  For  example,  if  we  take  to  be  the  free 
field  with  mass  0  =  m2 ,  and  d  >  2.  then  by  (3.5)  we  have  that 

=  ,->"+j-2(2t)  ||A||2+W 


+  0 


d\ 


after  the  obvious  change  of  variable.  If  this  is  to  be  equal  to  G'(.4).  so  that  by  Theorem  6.1 
will  be  self-similar,  we  require  a  =  \d-  1.  and  0  =  0.  The  condition  on  8  is  unpleasant, 
for  it  implies  that  only  the  massless  free  held  is  self-similar. 

From  the  point  of  view  of  Physics  this  is  undesirable,  since  it  is  well  known  that  all 
free  fields  are  renormalisable  (self-similar)  as  long  as  one  renormalises  mass  at  the  same 
time  as  renormalising  space.  That  is.  we  must  be  allowed  to  change  8  in  some  fashion. 
In  order  to  see  how  to  do  this,  it  is  most  natural  to  leave  the  spectral  setting,  and  think 
of  Gaussian  fields  as  arising  from  sums  of  Markov  paths.  Thus  for  the  remainder  of  this 
section,  we  treat  only  Gaussian  fields  whose  covariance  function  is  also  an  appropriate 
Green  function.  We  start  by  looking  at  the  structure  of  the  Markov  paths  themselves. 

6b.  Self-similarity  for  Markov  processes. 

Let  A ,  as  usual,  be  a  Markov  process  with  initial  measure  m(d.r)  =  Lebesgue  measure 
on  FJ.  We  call  A  selj-simlar  with  index  3  if  for  every  rj  >  0.  A'  and  the  process 

(6.4)  iX[t):=t,X[„-U) 

are  identical  in  distribution,  where  A  is  a  process  (on  another  probability  space)  which 
has  the  same  transition  probabilities  as  A  but  initial  measure  t^‘m[dx). 

Note  how  we  have  had  to  change  the  initial  distribution  in  order  to  preserve  the 
distribution  of  starting  points  after  scaling  A  by  the  factor  r>.  If  we  were,  for  example,  to 
start  all  our  processes  from  zero  no  such  condition  would  arise.  (See  Lamperti  (1972)  for 
a  theory  of  self-similar  Markov  processes  under  such  a  condition.) 

Lemma  6.1.  A  has  transition  density  p{t:  t.  y)  satisfying 

(6.5)  x,  y)  =  tfdp{t.  r]x.  i\y)  for  all  r/  >().<>  0. 
iff  it  is  self-similar  of  index  3. 


Proof. 

(i)  Sufficiency:  For  A  C  Bd. 


r'J’n(,)X[,r3t)  e  .4) 


=  /  yddx  /  dyp(y  st:x,y) 

Rd  JA/r) 

=  f  ndd*  [  dytidp(t\T}i,ny) 

Jrj  JA/n 

-  I  du  I  dvp{t\  u.  f) 

JlRd  JA 

=  Pm(X(t)  €  A). 

(ii)  Necessity:  By  self-similarity 


by  (6.5) 


P{n X{n  d0  e  A  I  1} X[y  stn)  =  x)  =  P{X(()  e  A  I  A'(<0)  =  *} 


i.e. 

Pr,^{t-to){^rl,V~lA)  =  Pt-to(x,A ), 

where  Pt(x,A)  is  the  obvious  transition  probability  function.  Setting  t  - 
transforming  the  last  equality  to  a  statement  about  densities  gives 

p{n~*s\ x,y)  =  t]dp(s\  r\x,  i}y) 


as  required. 

Lemma  6.2.  Let  g6  be  the  Green  function  of  X  (and  so  also  of  X).  If  X  is 
with  parameter  J  then  for  all  i)  >  0 

(6.6)  9$[*,u)  =  nd~*99v  Ur/t-iy). 


Proof. 


</”  '(nz.iv) 


-9 


-f 


p{t\qx,riy)exp(-tOn  g  )dt 


e  ** p(et}:)\  yx,  yy)ds 


y»CX/ 

=  n^~d  j  e~f>sp(s:  x,  y)ds  by  (6.5) 

=  1M/(m/), 


t0  =  s  and 


self-similar 


from  which  (6.6)  follows  trivially. 


As  an  example  is  probably  in  order  at  this  point,  the  reader  might  like  to  check  for 
himself  that  Brownian  motion  is  self-similar  with  .1  =  2.  Relation  (6.6)  is  most  easily 
verified  from  the  Fourier  transform  representation  (3.5)  of  the  Green  function. 

Note  that  whereas  Lemma  6.1  gives  a  necessary  and  sufficient  condition  for  self¬ 
similarity,  (6.6)  is.  as  it  stands,  only  necessary.  If.  however,  we  assume  (6.6)  for  all  6  >  0, 
then  it  is  easy  to  see  that  is  also  sufficient. 

6c.  A  prologue  to  renormalisation. 

The  above  shows  us  that  many  Markov  processes  are  self  similar.  It  thus  seems 
reasonable  that  if  we  sum  appropriate  functionals  of  these  processes  then  the  limiting 
sums  -  that  is,  the  Gaussian  fields  and  their  functionals  -  should  somehow  inherit  this  self¬ 
similarity.  There  is.  however,  one  more  parameter  in  this  picture  that  requires  normalising. 
Inherent  in  our  construction  of  Markov  functionals  was  an  exponential  damping  factor.  In 
constructing  p\  from  X.  the  time  change  is  compensated  for  via  a  change  of  initial  measure 
and  a  scaling  factor.  In  building  functionals  on  pX  that  will  be  distributional^  equivalent 
to  those  on  A',  we  shall  also  have  to  compensate  in  the  damping  factor  for  the  different 
time  rates. 

Another  way  to  think  of  this  is  to  recall,  as  we  noted  in  §2b,  that  the  exponential 
damping  is  essentially  equivalent  to  exponential  killing.  From  this  viewpoint  it  is  clear 
that  when  we  change  the  time  scale  to  obtain  pX  from  A',  we  must  also  change  the  rate 
at  which  we  kill  pX. 

We  need  one  more  result  before  we  can  continue,  that  relates  the  distributions  of 
functionals  of  X  to  those  of  pX.  As  one  would  expect,  for  the  right  sort  of  functionals, 
self-similarity  is  available. 

Theorem  6.2.  Take  X  a  right  Markov  process  with  symmetric  Green  function  ge .  Let 
k,  k  >  1,  and  define  p  =  -j(q.  t})  by 


(67) 


^(A)  =  rri+Jk-;(n'lA). 


Then,  if  X  is  self-similar  with  index  d, 

F*[Xi . Xk)=q~k(  9+'i),3F**-*v)[Xi . A'fc),  a  =  k(d  +  d)/2. 

Proof.  For  notational  convenience,  we  assume  that  •)  is  absolutely  continuous  with  density 
<?(x)  on  R'u.  The  general  case  is  handled,  as  usual,  by  passage  to  the  limit.  Then  (6.7) 
becomes,  with  tx  —  k[d  +  d)  I 2. 

(6.8)  £$(x)  =  r)~k{d+3)/2q(x/n) 

The  first  thing  we  must  check  is  that  ‘dq  6  M0,?  *’*;  i.e.  {^q,^q)en~ »  <  °°-  But,  by 


(2.10), 


<«W-*  =  J  anqMg6ti  a(^ymy)dy 

=  q~k{d+3]  J  q(xf n)g6r,">> {x,y)q{y fq)dxdy  by  (6.8) 
j  q{u)g6v  * (r)\i.  r)v)q{v)dud\ 

=  j  q{\i)g9(u,v)q(v)d\idv  by  (6.6) 

=  (?<?>*< 

which  is  finite  by  assumption.  Thus  F^g  ^  (  is  certainly  well-defined. 

The  next  thing  we  shall  need  to  know  is  what  happens  to  F8(Xi , . . . ,  A* )  when  we 
change  the  initial  distribution  of  the  A',  from  m(dx)  to  qdm(dx).  But  this  is  clear,  for  the 
characterisation  (2.11)  between  F 8  and  *;  gives  us  that 


(6.10) 


{F®(A'i . A*),  with  initial  distribution  m(dx)} 

={K  (A'i,  •  •  ■ ,  A'*)  with  initial  distribution  qdm(dx)}. 


This  is  all  the  background  we  need.  Now  note 


. A'fc)  =  f  e-8"  °q(XAti) . Xk(tk))dtl...dtk 

=  ,,-k{d+3)/2  f  eX8"  '•|*,9(^-iA'1(<1) . T}~lXk(tk))dtl  . .  .dtk 

J  R* 

t- 

iq-k{3~d  1/2  f  e~8r>  5ltl<7(f?-lA'l(/l) . rt-lXk(tk))dti...dtk, 

Jr" 
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where  the  A\  have  initial  measure  ri~dm(di),  and  we  have  applied  (6.10).  Map  t,  — 
in  the  last  integral  to  obtain 

. Xk)±i Vkl9+d)l%  J  . „~lX rk))dri . .  .drk 

LrjKl*+*)/*F»(Xl . A'jt) 

by  self-similarity  via  (6.4).  This  proves  the  theorem. 

We  are  now  in  a  position  to  return  to  the  study  of  Gaussian  fields. 

6d.  Renormalising  Gaussian  fields. 

The  basic  renormalisation  of  a  function  indexed  field  is  the  given  by  (6.1).  As 
we  noted  earlier,  if  9  is  called  self-similar.  We  wish  to  extend  the  notion  of  self¬ 

similarity.  by  starting  with  a  family  of  covariance  kernels  ge(x.y).9  >  0,  which  are  the 
Green  functions  of  a  Markov  process.  For  each  9,  V*  will  be  the  corresponding  Gaussian 
field.  We  call  the  family  {$*,  9  >  0}  renormalisable  with  renormalisation  parameters  (a,  r) 


(6.11) 


f. 


If  the  are  measure  indexed,  set 


(6.12) 


Wil  =  *127), 


where  “7  is  defined  at  (6.7),  in  order  to  define  the  renormalisation  parameters. 

This  notion  extends  that  of  self-similarity  in  two  directions.  Firstly,  the  two  fields 
need  not  be  defined  on  the  same  probability  space.  Secondly,  there  is  an  extra  parameter 
available  in  the  renormalisation  in  order  to  match  the  fields.  This  corresponds  to  mass 
renormalisation  in  quantum  field  theory  (c.f.  Glimm  and  Jaffe  (1982)).  and  is  completely 
missing  in  the  usual  theory  of  self-similarity.  We  have 

Theorem  8.8.  Let  .9  >  0.  be  a  Green  function  of  a  seif  similar  Markov  process,  with 
index  d,  and  let  {$9,0  >0}  be  a  family  of  Gaussian  Belds  with  covariance  kernels  {g°  .9  > 
()}.  Then  the  are  renormalisable  with  parameters  (( d  +  d)/2.  d). 


1 


Proof.  For  this  result  we  do  act  need  to  explicitly  consider  the  Markov  process  corre¬ 
sponding  to  g$ .  In  fact,  since  all  fields  are  mean  zero  and  Gaussian,  we  need  only  check 
that  covariances  match.  Thus,  take  f.h  €  S.t,  set  a  =  (d  +  d)j 2  and  note 


E{^  Uf)-yn'0(h)} 

=E{ 

=  J  I  *[x,v)tyiv)d*dy 

=  n~2Q+*~dJ  J  J{tltl)Qe{xh-yh)h{yli})dxdy  by  (6.2),  (6.6) 

=  //  f(u)gd(n.  v)h(r)dudv 
=  E{*9(f).**(h)}. 

Thus  the  covariances  do.  in  fact,  match,  and  the  theorem  is  proven. 

In  the  final  result  of  this  section  we  consider  the  renormalisability  of  additive  function¬ 
als  of  Gaussian  fields.  In  essence,  we  have  already  done  this,  for  we  know  via  the  specific 
results  of  §4  and  the  more  general  results  of  §8  that  all  additive  functionals  of  Gaussian 
fields  can  be  approximated  by  sums  of  additive  functionals  (c.f.  (4.5))  and  Theorem  6.2 
tells  us  when  the  latter  are  renormalisable. 

For  completeness,  however,  let  us  state  the  following  result  for  the  stationary  case,  for 
which  the  reader  does  not  yet  have  to  know  the  results  of  §8. 

Theorem  6.4.  Under  the  assumptions  of  Theorem  6.3,  take  a  7  6  Define  for 

T)  >  0.  “7  as  in  (6.7).  Let  ¥*(7)  and  ¥t,??  *(^7)  be  as  in  (4.4).  Then 

(6.13)  q  =  k(d  +  d)/2. 

Proof.  The  result  follows  trivially  from  Theorems  4.3  and  6.2.  We  need  only  check  that 
;)*,  G  But  since  -■  is  symmetric  with  inverse  Fourier  transform  of  the  right  kind, 

it  is  automatically  in  "7.  Thus  only  integrability  needs  to  be  checked.  This,  however,  is 
the  calculation  at  (6.9),  so  we  are  done. 

In  the  terminology  of  normalisation  parameters,  we  have 


7.  Moment  Proofs 

In  this  section  we  shall  give  a  proof  of  the  limit  theorems  of  Section  4  via  the  method 
of  moments.  For  reasons  we  shall  outline  below,  this  methodology'  works  perfectly  for 
Theorems  4.1  and  4.2,  which  treat  only  convergence  to  the  Gaussian  limit,  but  breaks 
down  in  Theorem  4.3  if  the  integral  (4.4)  defining  the  functional  is  more  than  two- 
dimensional.  (Or.  equivalently,  if  the  F 9  of  (4.5)  is  defined  on  more  than  two  X,  at  a 
time.) 

In  essence,  the  proofs  of  this  section  are  redundant,  for  more  powerful  techniques, 
which  work  in  more  general  situations,  will  be  used  in  the  following  section.  Nevertheless, 
we  include  them  for  two  reasons:  firstly,  they  are  conceptually  easier  and  considerably  more 
concrete  than  those  of  Section  8,  and,  secondly,  we  find  it  both  interesting  and  illuminating 
to  calculate  the  moments  of  the  limiting  distribution  per  se. 

To  make  life  a  little  easier,  we  note  the  obvious  fact  (e.g.  Reed  and  Simon  (1972),  p. 
51)  that  linear  combinations  of  functions  of  the  form 

(7-1)  fo(fi . **)=/(*,)•••/('*)  /  € 

are  dense  in  and  linear  combinations  of  measures  of  the  form 

(7.2)  f0[dxi,...,dxk)  =  i{dxi)...f(dxk) 

are  dense  in  MeSymm.  Consequently,  since  the  functionals  $*(•),  (•),  Vte(-)  and  ¥*(•)  are 

all  continuous  in  probability,  it  clearly  suffices  to  prove  the  Theorems  of  Section  4  for  f0 
and  70  of  the  form  (7.1)  and  (7.2)  respectively. 

In  this  section  we  shall  limit  our  discussion  even  further,  to  the  case  of  k  =  2,  because 
of  the  following  result. 

Lemma  7.1.  For  /  €  Sj  and  7  6  the  moments  of  o^f/)  -f  .l^h),  a.  li  e  R. 

determine  the  joint  distribution  of  Q9{f)  and  ’i'*  hi¬ 
proof.  By  Feller  (1979).  Problem  30.3,  for  example,  the  joint  distribution  is  determined 


by  the  mixed  moments  if.  for  some  a  >  0. 


_  £  f  I  y| n  \ 

(7.3)  limn_  v - '■  -  <  oo, 

with  A'  =  $s(/)  or  ) -  But  standard  results  on  multiple  Wiener-Ito  integrals  (e.g. 
Major  (1981),  p.  119)  tell  us  that  the  n-th  moment  of  a  fc-fold  integral  is  of  order 
exp[(fc/2)n log n  +  <9(«)].  Since  $6  and  correspond  to  k  =  1  and  k  =  2,  condition 
(7.3)  is  satisfied.  Note,  however,  that  for  k  >  2  this  condition  is  not  satisfied,  and,  indeed, 
moment  methods  do  not  work  in  this  situation. 

We  now  turn  to  calculating  the  mixed  moments  of  the  lemma,  for  which  we  shall 
need  certain  Feynman  diagrams.  The  diagrams  that  will  suffice  for  our  purposes  consist 

of  a  set  of  vertices ,  labelled  1 _ _  m,  and  legs  labelled  (i,j),  j  =  1  which  belong 

to  the  i-th  vertex.  We  shall  require  only  the  cases  fc,  =  1  or  2.  A  diagram,  denoted  by 
(ki,...,km)  is  formed  by  pairing  all  the  legs  in  such  a  way  that  the  legs  in  each  pair 
belong  to  different  vertices.  Figure  7.1  shows  the  two  unique  (up  to  permutations  of  legs 
belonging  to  a  given  vertex,  or  of  like  vertices  among  themselves)  diagrams  for  (1. 1,2.2). 
(Note,  that  in  the  notation  of  Major  (1981)  and  Dobrushin  (1979)  we  consider  only  what 
they  call  the  complete  diagrams.)  Note  that  in  general  there  are  K  =  +  . . .  +  km  legs, 

and  so  A'/2  pairs. 


i 

3  <  >  , .  >- 

Figure  7.1.  The  two  diagrams  for  (1.1. 2, 2). 

The  broken  lines  represent  bonds,  the  heavy  lines  represent  legs. 


We  can  now  turn  to  the  moments  of  +  i't1’.  Clearly,  it  will  suffice  to  find  an 
expression  for 

E{\*eif)\m\*eh)}n-rn} 

(7'4)  =E{[J  /(A)ZG.(rfA)]m.[£  J  J  7(Al)7(A2)ZG.(dA1)ZG#(dA2))"-"‘}. 

for  n  =  1,2 .  m  =  0,1 . n.  But  this  is  now  easy,  for  prepared  formulae  for  the 

expectation  of  moments  of  Wiener-Ito  integrals  already  exist.  To  present  them,  let  T(n.  rn) 
he  the  collection  of  all  n-vertex  Feynman  diagrams  for  which  the  first  m  vertices  have  only 
one  leg,  and  the  remaining  n  —  m  vertices  two  legs  each.  Since  there  are  m  +  2(n  —  m)  = 
2 n  -  m  legs,  and  so  n  -  ~m  bonds,  m  must  be  even. 

We  number  the  legs  by  an  index  a.  a  =  1, 2n  -  m,  and  denote  a  bond  by  a  pair 
(at.a2).  (Clearly,  not  all  2-tuples  are  possible.)  Let  A  be  the  set  of  a  for  which  the  a-th 
leg  comes  from  a  2-leg  vertex.  Define 

/|/(A)|2G*(dA)  a,d6.4c, 

/|7(A)|2G’s(dA)  c,3€A, 

J  7(A)r(A)G*(dA)  ae.4Je.4c 
or  or  €  Ac,  3  €  ,4. 

Note  that  in  the  third  case  the  various  symmetry  requirements  give  us  that  ras  is  also 
given  by 

(7.6)  j  i[\)f(-\)G»[d\)  =  j  7(-A)/(A)GVA). 

We  can  now  state 

Lemma  7.2.  For  f  €  Sd  and  7  €  *M^2mf7, ,  n  >  1.  m  =  0, 2, . . . ,  2|n/2), 

(7-T)  I  =  2"-"Er0,3i  . . .  ra,3, 

where  the  sum  is  taken  over  all  diagrams  in  r(n.m),  and  q  —  n  -  jm.  For  m  odd.  the 
mixed  moment  is  zero. 

Proof.  All  the  hard  work  has  already  been  done,  in  setting  up  the  notation.  The  Lemma 
is  a  straightforward  special  case  of  Corollary  5.4  of  Major  (1981). 

Before  we  can  finally  turn  to  the  proof  of  Theorem  4.3.  we  need  some  information  on 
the  moments  of  functionals  of  Markov  processes. 


Lemma  7.3.  Let  X  be  a  Markov  process  as  in  Section  2 ,  and  *,  i ,  *  2  6  M/>1.  Then 


■(F^(X)F^(X)}  =  D~l  I  M(djt)g*[z.y) i2(dy) 


In  particular, 


=  0  1  {-yi.  ^2>e- 


E{\Fe1(X)\i}=0~l(1,1)e. 


Lemma  7.4.  Let  X  be  a  Markov  process  as  in  Section  2,  and  7  €  M 8,1 .  Then ,  for  n  >  1. 


E{|F?|.V)n  =  ^  i./.w^i).. -1(^.1 

(7.9)  <  f  ■  j 

Both  of  these  lemmas  are  of  the  same  form  as  Theorem  5.2  of  Dynkin  (1984),  where 
similar  moments  are  calculated,  albeit  for  a  slightly  different  class  of  processes  and  func¬ 
tionals.  It  is  a  straightforward,  although  tedious,  exercise  to  follow  through  Dynkin’s 
calculations  to  obtain  the  lemmas. 

We  shall  also  require 

Lemma  7.5.  Let  X  be  a  Markov  process  as  in  Section  2,  and  72  €  i  such  that 
~l2(dxi,dz2)  =  7i(dx1)7i(dj,2)  with  7,  6  X8,1.  If  F83  and  F8(  are  the  corresponding 


additive  functionals,  then 


F?j(X,,A-j)=F'i(.Y,)F*(.Y2). 


Proof.  In  the  notation  of  (2.13)-(2.14),  we  have 

hl{X\,*2)l2{dXl,dX2)  =  72(<fci,cf*3) 

=  li(dxi)-n(dx2) 

=  M*i)M*2bi(d.ri)7i(rfjc2)- 

for  some  finite  7'.  Thus,  with  S  =  (<5i.<52). 

h2(/,./2)687j  ^(/|,  r2)  =  hi(x,)h,(/2)  J  e~ei'6'+s'l]  p6l  (r , .  y{  )p«,(/2.  y2hi  (dXi  )Va(djr2 ) 


Applying  the  above  equivalences  to  the  defining  form  for  the  F:  establishes  the  Lemma. 

We  can  now  finally  turn  the 
Proof  of  Theorem  4.3  for  the  case  k  =  2 

Firstly,  we  shall  need  to  add  one  more  assumption  to  those  surrounding  Theorem  4.3. 
viz.  the  integrals 


(7.10) 


J  .  .  .  J  g-  ( Xy .  4-2  )  •  •  •  ge  (  /«  - 1 ,  *n  b  (dll  )  .  .  . 


(d*n) 


are  finite  for  all  n  >  1.  Since  by  Lemma  7.4  these  integrals  bound  the  n-th  moment  of  the 
limit  variable  ®tf(7)  of  the  theorem,  their  finiteness  is  a  natural  ingredient  of  any  moment 
proof. 

We  commence  by  noting  that  by  Lemma  7.1  we  need  only  prove  convergence  of  mixed 
moments  to  obtain  the  convergence  in  distribution  of  the  Theorem.  By  (7.1),  (7.2)  and 
the  comments  there  we  can  restrict  ourselves  to  72  €  *A<symm  which  are  of  the  form  7x7, 
with  -  6 

To  make  life  notationally  easier,  let  us  write  the  field  $*(/)  of  (4.1)  in  the  form  $^(7, ) 
of  (4.3),  where  7!  =  7t(/)  is  defined  by  7 ,(.4)  =  JAf(x)dx ,  and  we  have  77  €  At*1. 
Consequently,  we  are  interested  in  studying  the  mixed  moments  of  ($^(71),  ♦^(72)).  with 
7 1  G  and  72  €  and  of  product  form.  We  have,  for  n  >  1,  m  <  n.  A  >  0, 

£{[^(7j)r[^(72)r-> 

/v(a)  _ 

(7iJl)  ‘=‘  i<;<fc<.v(M 

=2m'n(0/\)n-l'*m  f;  ^^E,£{^(i.j,k)}.£{^(i,j,k)}, 

;V=2 

where  i  =  (i, . im),  j  =  (;, . jn-m)-  k  =  (kj . k„_m), 

^'T(i.j.k)  =<r,,  ...(rlm<ru  . . .  <r]n  _m<rkx  ...<Tkf)_m, 

rrtl-l  k)  =  K (X, )  .  K  (XJK(.X„ .  X„,  | . . .  J 

=  F> ,  ( .Y, ,  |  . .  J*  ( X.  „ )  F *  |  .Y, , ) . . .  J*  ( .Y,„  _  |  F*  ( X„ , ) . . .  C?  ( A* . .  „ ) . 


(the  last  equality  following  from  Lemma  7.5),  and  the  sum  E'  in  (7.11)  is  over  all  i.j.  k 

with  components  in  (1 . A’}  such  that  ja  7 -  ka,  ct  =  1 . ti  —  m. 

Consider  the  expectation  of  the  product  7rCT(i. j,  k).  This  is  clearly  zero,  unless  each 
(ra  appears  an  even  number  of  times  among  the  2 n  -  m  <r’s.  Thus,  if  m  is  odd,  the  mixed 
moment  is  zero,  as  is,  by  Lemma  7.2,  is  the  mixed  moment  of  and  4,e.  Henceforth, 
therefore,  we  assume  m  is  even. 

We  now  consider  those  configurations  (i.j,  k)  composed  only  of  q  =  n  -  m/2  different 
indices,  (ai, . ..  ,a„),  each  appearing  precisely  twice.  Then,  clearly. 

9 

(7-12)  »<T(U,k)  =  n<sL 

i/=i 

Consider  the  corresponding  fff(i.j,k),  appearing  in  (7.12).  To  each  <t\v  in  (7.12)  we  can 
associate  a  term  rQi,  coming  from  the  expectation  of  the  product  of  two  functionals  of  XQi.. 
according  to  the  following  rules:  (c.f.  Lemma  7.3).  If  a  =  iu  =  jb.  or  o  =  ia  =  kb .  then 

(713)  r„  =£{J^,(.\-„JF*(.V0  J)  7)„. 

If  a  =  ia  =  ift,  then 

(7  U)  'a  =  £{J*  (-V„Jf1,(X,,)}  = 

If  a  =  ja  =5  jb,  or  a  =  ja  =  kb,  or  a  =  ka  =  kb,  then 


(7-15)  rQ  =£{Fe7(XaJF?(A’aJ}  =  (7.  7),. 

Consequently,  still  in  the  case  of  precisely  q  different  indices,  the  expectations  of  the 
irF  terms  in  (7.11)  depend  only  on  the  patterns  formed  among  the  indices,  and  not  on 
the  indices  themselves.  Now  note  that  we  could  choose  (oi,...,a?)  from  (1,...,.V)  in 
(Ar!)/((A’  -  g)!}  different  ways,  so  that  we  have 


(7.16) 


=  2m~n  2 ^ 


Nl  (A r-q)'. 


- —A  q T.'6qE{nF (i,  j,  k) }  +  other  terms 


_  2m~nT,'0'1  E {irF (i,y  k)}  +  other  terms, 

where  the  sum  is  over  all  the  (i,j,k)  of  (7.11)  of  the  paired  form  discussed  above.  The 

"other  terms”  come  from  choosing  q1  <  q  different  indices  (qj . a^).  so  that  some  of 

the  (T,  appear  2 p  times,  p  >  1.  These  terms  are  asymptotically  negligible,  since  under 
(7. 10)  boundedness  of  the  moments  of  and  F*  (Lemma  7.4)  guarantee  that  the  (q  -  1) 


“other  terms”  are  of  order 


A=7'  *  '  '*  H  ’’ 


Thus  it  remains  for  us  to  consider  the  A-independent  expression 


(7.17)  2m~nF'6Q {rrF ({,},  k)}. 

This  is  clearly  expressible  in  terms  of  the  rQ  of  (7.13)-(7.15).  Now  recall  the  initial  defini¬ 
tion  of  7 1  in  terms  of  /.  to  note  that  if  Ge  is  the  Fourier  transform  of  g e  then: 


<7„1.).=  /  |/(A||»G*(dA). 
(7.7).  =  j  l7(A)|2G»(<iA), 

(7.7,).  =  J  7|A)/’(A)GVa). 


That  is,  the  ra's  of  ( 7, 13)— ( 7.15)  correspond  to  the  rag' s  of  (7.5),  modulo  a  factor  of  6. 
Furthermore,  we  can  write  the  sum  in  (7.17)  as  a  sum  over  Feynman  diagrams  r(n,  m),  in 
which  the  m  single-leg  vertices  correspond  to  indices  in  i.  and  the  n-m  double-leg  vertices 
correspond  to  indices  in  j  and  k.  Thus,  since  the  factors  of  8  cancel  nicely,  we  have  that 
(7.17)  is  precisely  (7.7).  That  is,  the  mixed  moments  of  $*(/)  and  4^(72)  converge  to 
those  of  $d(/)  and  4^  (72).  This  completes  the  proof. 


A  remark  on  k  >  2.  It  is  possible  to  show  that  even  if  k  >  2  the  moment  sequence  of  the 
approximating  sums  converge  to  those  of  the  Gaussian  limit.  The  only  essential  difference 
in  the  proof  is  that  more  complicated  Feynman  diagrams  are  required.  Wolpert  (1978b) 
has  actually  done  this  for  the  special  case  of  A'  Brownian  motion  in  and  the  additive 
functional  of  the  limiting  field  a  Wick  power.  Note,  however,  that  for  k  >  2  the  moments 
of  Gaussian  functionals  do  not  determine  their  distributions,  and  so  no  convergence  in 
distribution  result  can  be  obtained  this  way. 

We  conclude  this  section  with  some  brief  comments  on  the 
Proofs  of  Theorems  4.1  and  4.2.  It  is  clear  that  both  of  these  results  are  virtually 
standard  central  limit  theorems,  so  that  all  we  need  check  is  that  the  summands  satisfy 
appropriate  moment  conditions  and  the  mean  and  variance  of  the  limit  are  correct.  Both 
of  these,  however,  follow  from  Lemma  7.3,  and  so  the  proofs  are  trivial. 


8.  The  General  Limit  Theorem 

In  this,  essentially  final,  section  we  shall  state  and  prove  a  more  general  result  than 
those  of  Sections  4  and  7.  The  motivation  for.  and  intuitive  appeal  of  this  result  is  as 
for  the  earlier  ones.  Having  seen  how  the  arguments  of  Sections  5  and  6  apply  to  the 
Theorems  of  Section  4,  the  reader  who  has  got  this  far  will  see  how  they  also  apply  in  the 
most  general  situation. 

In  order  to  set  up  this  general  result,  we  shall  need  a  theory  of  multiple  Wiener  inte¬ 
grals  far  more  general  than  that  of  Section  3.  We  thus  recall  such  a  theory,  due  to  Dynkin 
and  Mandelbaum  (1983).  that  has  its  roots  in  Neveu  (1968).  We  follow  Mandelbaum  and 
Taqqu  (1984)  most  closely. 

8a.  Wiener  integrals  on  an  arbitrary  space. 

Let  (X,3,is)  be  an  arbitrary  measure  space,  and  define  a  Gaussian  process  /j  on 
C2(v)  with  zero  mean  and  covariance 

(8.1)  £{/i(/i)/i(/2)}  =  W/i/2):=  f  fiix)ft{x)v{dx). 

If  we  consider  the  subfamily 

{&'(£)  =  /i(lfl),  B  €8,v(B)  <oc} 

as  a  Wiener  process  on  X,  then  It  can  be  written  symbolically  as 

/.(/)  =  J  /(*)»•(<**). 

Now  write  Mk  =  £]yrnm{vx- .  ,xy)  for  the  space  of  symmetric  functions  hk(z  i,  —  -t*) 
for  which 

uk(hl)  =  J  hl(zl,...,xk)i/(dzl)...i/(dzk)  <  oc. 

The  multiple  Wiener  integral  of  order  k  is  a  linear  mapping  Ik  from  Mk  into  the  space 
of  random  variables  which  are  functionals  of  the  Gaussian  family  Ii(f).  The  mapping  is 
defined  uniquely  by  conditions  .4  and  B  following: 


(A)  For  functions  of  the  form 


(8.2)  h[{xy . xk)  =  /(/,).  ../(**!.  f€£*M 

we  have 

(8.31  /*(*£>  = 

where  is  the  Hermite  polynomial  of  degree  k  with  leading  coefficent  1. 
(B)  For  hk  €  Mk, 

E\Il(hk)\=k\vk[hl). 


Again,  symbolically,  it  makes  sense  to  write 


Ik(hk) 


-L  L 


hk{x\ . xk)W(dxl) ..  .W(dxk). 


It  is  important  to  note  for  latex  use  that  the  linear  combinations  of  functions  satisfying 
(8.2)  are  dense  in  Uk)  and  so  their  integrals  are  dense  in  the  space  of  all  integrals. 

8b.  A  general  limit  theorem. 

The  following  theorem  follows  from  Theorem  2  and  the  discussion  of  Section  1.5  of 
Dynkin  and  Mandelbaum  (1983). 


Theorem  8.1.  Let  A’l.A'a....  be  independent  and  identically  distributed  random  vari¬ 
ables  taking  values  in  ( X,  B ),  and  with  distribution  v.  For  A  >  0  let  N \  be  a  Poisson 
variable  with  mean  A,  independent  of  the  A,.  Let  be  an  i.i.d.  sequence  of 

variables  with  P{a,  =  +1}  =  P{<r,  =  -1}  =  and  independent  of  both  the  N\  and  A*,. 
Finally ,  for  k  =  1,2, ... ,  take  hk  6  Mk  and  let  Ik(hk)  be  its  multiple  Wiener  integral  as  in 
the  previous  subsection.  Then  the  random  variables 

OG 

(8.4)  Zx{hl'ht...):=£\-k'9  £>„..  .<7,tMA'„ . A',J. 

k-l 

with  the  inner  sum  over  1  <  »i  <  »2  <■•■<»*  <  .VA,  converge  in  distribution ,  as  A  — •  oo, 
to 

Z(hlth2...):=  T±Ik[hk). 


(8.5) 


As  noted,  an  elegant  proof  of  this  powerful  result  can  be  found  in  Dynkin  and  Man- 
delbaum  (1983).  In  Mandelbaum  and  Taqqu  (1984)  the  result  is  extended  to  an  invariance 
principle. 

Now  let  us  consider  a  special  case  of  Theorem  8.1.  Let  X  be  a  right  Markov  process 
with  symmetric  Green  function  g6.  Let  _T  =  (R“)r  +  be  the  path  space  of  A',  and  u  the 
measure  the  process  induces  on  X.  For  k  >  1,  and  7*  £  Mk'0 ,  let  F%t  be  the  corresponding 
functional  (as  in  §2.b)  of  k  independent  copies  of  X.  Then,  in  the  above  notation,  F^k  €  Uk-, 
and  so  the  multiple  Wiener  integral  Ik(F^)  is  well  defined.  Since  there  is  a  one-one 
correspondence  between  functionals  i*(F®  )  and  the  pair  (7*,#),  we  save  on  notation  by 
setting 

(8.6)  lkhkie)  =  Ik(Fl). 


Let  X\  and  <rx-  <r2. . . ..  be  as  in  Theorem  8.1,  and  Ax,  X2, . . .,  i.i.d.  copies  of  X.  We  have, 
as  an  immediate  consequence  of  Theorem  8.1, 

Theorem  8.2.  For  7 k  G  k 


£(»/A)‘/2I . A-,. I 

k=  1 

converges  in  distribution,  as  A  —  00,  to 


00  flk/2 
k=l 


This  result  completely  incorporates  the  results  of  Section  4.  and  is,  in  fact,  far  more 
general,  since  far  fewer  restrictions  are  placed  on  the  7*.  The  restrictions  there,  that 
7 k  €  *A<symm*  were  placed  primarily  to  permit  easy  identification  of  the  limit  as  a  simple 
Wiener-Ito  integral.  In  the  more  general  result  above,  it  is  not  quite  as  clear  what  the 


IkhkJ)  are. 

To  relate  Theorem  8.2  with  what  has  gone  before,  let  us  use  it  to  obtain  a  new 
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Proof  of  Theorem  4.S. 


We  have  to  prove  convergence  of  the  pair  ($^(/).  ^(7))  of  (4.1)  and  (4.5)  to  the  pair 
{$*(/).  ¥*’(7))  of  Section  3  and  (4.4).  The  restrictions  of  Theorem  4.3  placed  /  in  S(</^). 
and  7  in  *Mds-ykmm. 

Firstly,  to  make  comparisons  easier,  affix  a  suffix  k  to  7:  i.e.  7*.  =  7.  Now  use 
7i  =  7i(/)  to  denote  the  measure  in  M9’1  with  density  /.  Then,  comparing  (4.1)  and 
(4.5)  with  the  definition  of  an  additive  functional  in  Section  2,  it  is  immediate  that  we 
have 

*'(/)  3  *1(7, (/I). 

Consequently,  by  Theorem  8.2  and  the  Cramer-Wold  device,  we  have  that  ($^(/),  ¥^(7*)} 
converge  in  distribution  to  some  pair  {01^2/i(7i),  *7£T-/fc(7fc,0)).  We  need  only  check  that 
this  is,  distributionally,  the  right  pair. 

We  know  that  /,(•)  is  Gaussian,  with  zero  mean.  Its  covariance  function  is  given,  for 
Vi  ,1*2  €  M9A,  by 

E{IMh  M}  =  e{f9if9j 

=  &  lf  /  Vi(dz)ge(x,y)v2(dy). 

(c.f.  Lemma  7.3).  Since  ~h  has  density  /,  the  last  integral  leads  to  var(/1(71))  = 
e~l  /  /  f[*)sl9(*,v)f[v)dzdy,  and  so  Ji(7i)=0_1  /2 Q9[f).  To  complete  the  proof,  we  con¬ 
sider  7 k  €  *M9Symm  of  the  form 


7 k{dxi . dxk )  =  7(d/i)...7(d-tfc),  7  € 

The  linear  combinations  of  such  7*  are  dense  in  and  so  it  will  suffice  to  consider 

them.  We  have 


(8.7) 


hUJ)  =  h{F9k) 


by  (8.3).  But.  since  the  equivalence  between  Ix  and  has  already  been  established,  and 
'  €  we  have 

(8-8)  J  7(A)ZG*(dA), 

via  the  Wiener-Ito  integrals  of  Section  3.  Now  apply  Ito’s  lemma  (e.g.  Major  (1981),  p. 
30)  to  (8.7)  via  (8.8)  to  obtain 

Hk(  /•••  / 

V^/eKFJP)/  J  1 

Thus 

0k/i  t  1  f  f 

—h  hk.«)=p  J  ■I  7(A1)...7(At|2G.(<<A1)...Zo.(<iA0 

=  *'00 

by  (4.4).  Furthermore,  the  above  argument  also  establishes  that  the  joint  distribution  of 
the  limit  variables  is  the  correct  one,  so  we  are  done. 


9.  Concluding  Comments. 

Some  comments  on  what  we  have  not  done: 

1.  The  connection  with  Physics.  Throughout  the  paper  we  have  made  numerous 
references  to  Euclidean  quantum  field  theory,  without  ever  being  really  precise.  This  was 
to  avoid  trebling  (at  least)  the  length  of  an  already  overlong  paper.  The  reader  now 
interested  enough  to  find  out  more  about  this  could  start  with  Glimm  and  Jaffe  (1982). 

2.  The  discrete  case.  We  have  not  considered  lattice  indexed  Gaussian  fields.  It  is  clear 
that  this  could  be  done  by  representing  them  as  the  sum  of  random  walks  on  the  lattice, 
rather  than  the  more  general  Markov  processes  we  have  considered.  Some  of  these  results 
could  be  obtained  from  ours  by  using  the  “discretization"  trick  of  Dobrushin  (1979).  This, 
however,  would  be  like  using  a  sledgehammer  to  kill  the  proverbial  ant,  since  a  direct 
approach  would  be  very  simple.  For  an  idea  of  how  this  might  go,  see  Williams  (1973), 
Spitzer  (1974).  or  Dynkin  (1983). 

3.  No  transition  density.  Not  every  Markov  process  has  a  transition  density,  nor  does 
every  Gaussian  field  have  a  covariance  functional  that  can  be  written  via  a  kernel.  In 
the  more  general  case,  Dynkin's  theory  of  additive  functionals  (§2)  based  on  the  Green 
function  and  symmetry  could  be  replaced  by  a  theory  based  on  the  resolvent  operator  and 
duality  (c.f.  Dynkin  (1982)  and  Getoor  and  Sharpe  (1985)),  and  so  a  more  general  class 
of  Gaussian  fields  considered.  Since  no  new  insight  would  be  forthcoming,  while  a  lot  of 
very  difficult  Mathematics  would,  we  shall  leave  this  to  our  betters. 
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